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Abstract 

We aim to construct a general framework for portfolio management 
in continuous time, encompassing both stocks and bonds. In these 
lecture notes we give an overview of the state of the art of optimal bond 
portfolios and we re- visit main results and mathematical constructions 
introduced in our previous publications (Ann. Appl. Probab. 15, 
1260-1305 (2005) and Fin. Stoch. 9, 429-452 (2005)). 

A solution of the optimal bond portfolio problem is given for gen- 
eral utility functions and volatility operator processes, provided that 
the market price of risk process has certain Malliavin differentiability 
properties or is finite dimensional. 

The text is essentially self-contained. 
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1 Motivation 

The literature on portfolio management starts with the Markowitz portfolio 
and the CAPM ([19], [20], [33]). It is a one-period model, where the infor- 
mation on assets is minimal. Every asset is characterized by two numbers, 
its expected return and its covariance with respect to the market portfolio. 
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With such poor information, one cannot hope to distinguish between stocks 
and bonds, and indeed part of the beauty of the CAPM hes in its generahty: 
it apphes to any type of financial assets. 

On the other hand, as soon as one tries to make use of all the information 
available on assets, important differences appear between stocks and bonds. 
Bonds mature, that is they are eventually converted into cash, whereas stocks 
do not. The price of bonds depends on interest rates, and the price of stocks, 
at least in the academic literature, does not. The bond market is notoriously 
incomplete, much more so than the stock market, as is observed in practice. 
As a result, the classical results on portfolio management, such as Merton's 
([21], [22]), concern stock portfolios. This paper and the papers [10] and [M] 
were born from a desire to extend them to bond portfolios. 

More generally, we aim to construct a general framework for portfolio 
management in continuous time, encompassing both stocks and bonds. 

The first difficulty to overcome (and, in our opinion, the main financial 
one) is the fact that such a theory should encompass two very different kinds 
of financial assets: bonds, which have a finite life, and stocks, which are per- 
manent. We do it by introducing a new type of financial asset, the rollovers. 
A rollover of time to maturity x is a bank deposit and which can be cashed at 
any time, with accrued interest, provided notice be given time x in advance. 
Roll-overs have constant time to maturity (as opposed to zero-coupon bonds, 
for instance), and are similar to stocks, in the sense that their main char- 
acteristics do not change with time. By decomposing bonds into rollovers, 
instead of decomposing them into zero-coupons, we can hope to incorporate 
bonds and stocks into a unified theory of portfolio management. Rollovers 
were considered in [32] under the name "rolling-horizon bond" . 

This implies that the time to maturity x, rather than the maturity date T, 
becomes the relevant characteristic of bonds. Thus, we shall describe bonds 
using a moving maturity-time frame, where at time t, the origin is the time 
to maturity x = 0, corresponding to the maturity date T = t. As we shall 
see very soon, there will be a mathematical price to pay for that. 

At any time t, denote by Pt{x) the price of a unit zero-coupon with time 
to maturity x. The function x ^ Pt{x) will be called the zero-coupon (price) 
curve at time t; note that the actual time when that zero-coupon matures 
is T = t + X, and that T, is fixed while x changes with t. The zero-coupon 
curve pt will be understood to move randomly, and the second difficulty we 
face is to describe its motion in some reasonable way. One solution is to 
decide that pt belongs to a fixed family of curves, depending on finitely many 
parameters, so that 

Ptix) = /(t,x;ri, ...,rd) 
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and the random motion of pt is the image of a random motion of the r^, 
which could, as in spot-rates models, be modelled, for instance by diffusions. 
This is the parametric approach, which exhibits the classical difficulty of all 
parametric approaches, namely that there is no theoretical reason why the 
Pt should be written in that way, so that the choice of the function / has 
to be dictated by observational fit. One then has to strike the right balance 
between two evils: if the number of parameters is too small, the model will 
be unrealistic, and if it his higher, it becomes very difficult to calibrate. 

We will operate in a non-parametric framework: we will make no assump- 
tion on Pt, beyond some very rough ones, regarding smoothness and behavior 
at infinity, nothing that would much constrain their shape. Mathematically 
speaking, we will let the curve pt move freely in a linear space E, which will 
typically be an infinite-dimensional Banach space, of functions from [0,oo[ 
to R. 

In order to reflect adequately known financial facts, the correct definition 
of E must incorporate some basic constraints: 

1. At any time t, the zero-coupon prices pt{x) must depend continuously 
on the time to maturity x. In order for forward interest rates to be 
well-defined, they must also have some degree of differentiability with 
respect to x. So E must consist of continuous curves with some degree 
of differentiability. 

2. The degree of differentiability of functions in E will determine which 
basic interest rates derivatives can be modelled. If pt is continuous, 
for instance, then we can introduce bonds. The price of a unit zero- 
coupon bond with time to maturity x is Pt{x)\ the bond itself, i.e. the 
value of a portfolio including exactly one bond is represented by the 
linear form pt i— >■ Pt{x). Mathematically speaking, this is just the Dirac 
mass Sx at x. Now other derivatives such as Call's and Put's on zero- 
coupon bonds can be introduced, since the pay-off for each of them is a 
continuous function of the zero-coupon bond price Pt{x), with a given 
time to maturity x. If pt is continuously differentiable, then the forward 
interest rate with time to maturity x, —^pt{x)/pt{x) is well-defined, 
and further contingent claims can be defined, such as caps, floors and 
swaps. 

3. The curve pt will be understood to move randomly in E, the random- 
ness being driven by a Brownian motion. We will therefore need to 
define Brownian motions in the infinite-dimensional space E, which for 
all practical purposes will require to be a Hilbert space. 
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4. The accepted standard in mathematical modeUing of zero-coupon prices 
(the Heath- Jarrow-Morton model, henceforth HJM) is to decide that 
the real- valued process t i— > pt{T — t), the price at time t of a unit 
zero-coupon maturing at a given time T, is an Ito process satisfying an 
stochastic ordinary differential equation (SODE). As is well-known, for 
fixed X, the real- valued process 1 1— > Pt{x), which is also an Ito process, 
then no longer satisfies an SODE. Indeed, if f(t,T) = pt(T — t), then 
we have pt{x) = f{t,t + x) so that for fixed x : 

dmix) = [dJ{t,T) + ^l^dT]T=t+. = [dJ{t,T)]T=t+. + ^^dt. 

(1) 

Here the right-hand side (r.h.s) depends, not only on pt{x), but also on 
its partial derivative with respect to x. So, equation ([1]) for p is a SPDE, 
stochastic partial differential equation, where the first term on the r.h.s 
depends only on the un-known pt{x), since f{-,T) satisfies an SODE. 
This is the well-known difficulty of the Musiela parametrization (see 
|25j). and the space E shall permit a simple mathematical formulation 
of the SPDE ([ID. 

5. At any time t, the zero-coupon prices Pt{x) should go to zero as the 
time to maturity x goes to oo. To include also the trivial case, where all 
interest rates vanish, and also cases where the forward rates converges 
rapidly to zero as x oo, we only require that \im pt{x) exists as 
X oo. N.B. We will chose E such that the elements f & E satisfying 
\im.x^oo /(^) = form a closed sub-space of E, in order to cover easily 
the case where pt{x) 0. 

Formula ([I]) is really an infinite family of coupled equations, one for each 
X > 0, describing the motion of the random variable pt{x), which we write 

dpt (x) 

dpt{x) = pt{x)mt{x)dt + pt{x)at{x)dWt H g^^t, (2) 

where for the moment W is thought of as being a high dimensional Brownian 
motion. Let us rewrite it as a single stochastic evolution equation for the 
motion of the random curve pt in E, SODE in E : 

dpt = Ptmtdt + Pt(rtdWt + {dpt)dt (3) 

where d is the differentiation operator with respect to time to maturity, i.e. 
it is defined by {du)(x) = for differentiable u ^ E. Since the left-hand 

side "belongs" to E, so must the right-hand side, and then dpt must belong to 
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E. There are ways to achieve that. One is to choose a framework where the 
operator d is continuous over all of E. Then so is its n-th iterate d'^, so that 
the space E must consist of functions which have infinitely many derivatives. 
Unfortunately, the natural topology of such spaces cannot be defined by a 
single norm, except for very particular cases, and the mathematics become 
more demanding. A second more standard way to proceed is to consider d 
as an unbounded operator in a Hilbert space E, so that d is defined only on 
a subspace T>{d) C E, called the domain of the operator. One would then 
hope to define the solution of equation ([3]) in such a way that, if the initial 
condition po lies in T>{d), then pt remains in 1^(19) for every t, so that t ^ pt 
is a trajectory in T>[d). In other words, if po{x) is differentiable with respect 
to X, so should the functions x t— > pt{x) be for all t > 0. 

To summarize, the introduction of rollovers and a moving frame forces 
us to complicate the equations for price dynamics, by incorporating an ad- 
ditional term, dpt- To be able to solve the relevant equations, we have to 
treat d as an unbounded operator in Hilbert space. The definition of the 
relevant Hilbert space has to incorporate basic properties which we expect 
of zero-coupon curves. 

This suits our purpose well, for it enables us to work in a non-parametric 
framework, where no particular shape is assigned to the the zero-coupon 
curves. On the other hand, we then have to use the theory of Brownian 
motion in infinite-dimensional Hilbert spaces and the corresponding stochas- 
tic integrals, which creates some additional difficulties. We do not limit 
the number of sources of noise, indeed in our paper there can be infinitely 
many. This is natural, since the already mentioned experimental fact, that 
even using a large number of bonds, not all interest rate derivatives can be 
hedged. The third difficulty to overcome, is the mathematically significant 
fact that such a market can not be complete in the usual sense, i.e. every 
(sufficiently integrable) contingent claim being hedgeable. This has impor- 
tant implications for the solution of the portfolio optimization problem. The 
now classical two-step solution, so successfully applied to the case of a finite 
number of stocks (cf. [17], [28]), consisting of first determining the optimal 
final wealth by duality methods and then determining a hedging portfolio, 
does not (yet at least) apply to the general infinite dimensional bond mar- 
kets. In this paper (see [TO] and [M]) we give, within the considered general 
Ito process model, the optimal final wealth for every case it exists (Propo- 
sition HS]). The existence of an optimal portfolio, is then established by the 
construction of a hedging portfolio for two cases : The first is for determin- 
istic i?-valued drift m and volatility operator cr, where we give a necessary 
and sufficient condition for the existence of an optimal portfolio. Here there 
can exist several equivalent martingale measures (e.m.m.), so the market can 
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clearly be incomplete in every sense of the word. The second is for certain 
stochastic m and a, for which there is a unique market price of risk process 
7. There is then a unique e.m.m. Q. Now, certain integrability conditions 
on the £^-valued Malliavin derivative of the Radon-Nikodym density dQ/dP 
leads to the construction of a hedging portfolio. 

We have tempted to make these notes self-contained, with exception of 
the general hedging result in Theorem [HHl The notes first recall some basic 
facts concerning linear operators and semi-groups in Hilbert spaces, Sobolev 
spaces and stochastic integration in Hilbert spaces. The theory of bond port- 
folios and hedging of interest rate derivatives are then introduced. Once this 
theory is explained, the paper proceeds to a short solution of the optimiza- 
tion problem, leading to the results of [TU] and In particular, under the 
assumption that the market prices of risk are deterministic, some explicit 
formulas are given, very similar in spirit to those who are known in the case 
of stock portfolios, and a mutual fund theorem is formulated. We conclude 
by stating an alternative formulation, of the optimization problem, within a 
Hamilton- Jacobi-Bellman approach. 



2 Mathematical preliminaries 
2.1 Hilbert spaces and bounded maps 

We shall be working with separable infinite-dimensional real Hilbert spaces. 
Let ii^ be a Hilbert space with scalar product { ■, )e and norm || ||^, simply 
denoted ( , ) and || || if no risk for confusion. The topology and convergence 
in E is w.r.t. this norm, if not otherwise stated, i.e. the strong topology and 
convergence. By definition E is, separable if it has a countable dense subset. 
One shows easily that E is separable iff it has a countable orthonormal basis 
e„, n G N, i.e. (cj, Cj) = for i 7^ j and ||ej|| = 1, so that every x & E can be 
written: 

00 

X ^ ^ i^X, Gri) 
n=0 

where the right-hand side converges in E. Since the e„ are orthonormal, we 
have Parseval's equality: 

00 

n=0 

A typical separable Hilbert space is i"^, which is the space of all real sequences 
an,n G N, such that X]|a„|^ < 00. The scalar product in i"^ is given by 
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{a,b) = 'Y^anhn- In fact, every infinite dimensional separable Hilbert space 
E is isomorphic to The map 

X t-^ a„ = (x, en)E, G N, (4) 

of E into £^ is a linear bijection and it preserves norms on both sides. 

A linear map L : £"1 — *• £"2 is continuous if and only if it is hounded, that 
is if there exists a constant c such that ||La;||^2 < c||a;||^^ for every x G Ei. 
The (operator) norm of L is then defined to be the infinimum of all such c: 

||L|| = inf {c I \\Lx\\j^^ < c||x||^^ Vx} . 

For example, the linear map in (jll) of £ onto £^ as well as its inverse has norm 
1. The linear space of all continuous linear maps from Ei to £2, £(£'1, £2), 
is a Banach space when given this norm. One writes L{E) as a shorthand 
for L{E,E). Linear maps are also called linear operators or just operators. 
A bounded operator on £ is a bounded linear map from £ into itself. The 
dual space £' of £, i.e. the space of all linear continuous functionals on £, 
is given by £' = L(£,]R). By the F. Riesz representation theorem, 

£ G £' iff 3/ G £ such that £(x) = (/, x) Vx G £. (5) 

Also ||£||£;' = so E' and £ are isomorphic. In this paper we will often 

use, in the context of Sobolev spaces, other representations of the dual £'. 

By duality, every operator in £(£i,£2) corresponds to an operator in 
£(£2, E[). Using the representation of the dual space given by ([5]), the adjoint 
operator A* oi A ^ L{Ei, £2) is defined by A*y = y*, where for y E E2 the 
element y* G £1 is defined by 

{y\ x)e^ = {y, Ax)e2 Va; G £1. (6) 

This defines an operator A* G £(£2,£i). One easily checks that {A*)* = A 
and 11^*11 = ll^ll . Let us consider a simple example, which will be relevant 
in the sequel of this paper: 

Example 1 (Left-translation in L^) 

i) Let £ = L^(]R) and let a be a given real number. Define the operator A 
on £ by iAf){x) = f{x + a). Then \\A\\ = 1 and iA*f){x) = f{x - a). We 
note that A has a bounded inverse A~^ given by {A~'^f){x) = f{x — a), so 
AA* = A* A = I, where I is the identity operator. 

ii) Let £ = £^([0,oo[) and let a > be a given real number. Define the 
operator A on £ by {Af){x) = f{x + a). Here we find that \\A\\ = 1, that 
a.e. {A*f){x) = OifO<x<a and that {A*f){x) = f{x — a) if a < x. In this 
case A* is one-to-one and AA* = I. But A* A is the orthogonal projection on 
the (non-trivial) closed subspace of £ of functions with support in [a, oo[. 

So A* A I. 
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An operator S G L{Ei,E2) is called unitary if SS* = S*S = I. This is the 
case of A in (i) of Example [H An operator S G L{Ei, E2) is called isometric 
if S*S = I. This is the case of A* in (ii) of Example [B 

We will be interested in a particular class of bounded operators on E. We 
begin with an easy result 

Lemma 2 Suppose L G L{Ei,E2) and that we have: 

00 



|2 ^ 
< 00 



n=0 



for an orthonormal basis Cn, G N in Ei. Let /„, n G N be another orthonor- 
mal basis. Then: 

00 00 

|2 



n=0 n=0 

Definition 3 An operator L on Ei into E2 is Hilbert- Schmidt if Yl'^=o ll-^^nll^ < 
00 for some orthonormal basis en-,n G in Ei. Its Hilhert- Schmidt norm is 
defined to be: 

> 1/2 

\^\\-HS = I E 11"^*^'^ 
\n=0 

(iocs not depend on the choice of the orthonormal basis en,n G N, in E. 
The linear space of Hilbert- Schmidt operators from Ei into E2 is denoted 
T-CS{Ei, E2). 

Hilbert- Schmidt operators are bounded (in fact, ||L|| < and even 

compact: they map bounded subsets of Ei into relatively compact subsets 
of E2. In other words, if L is Hilbert-Schmidt and (a;„)„gj^ is a bounded se- 
quence, then one can extract from (La:„)^gp^ a norm-convergent subsequence. 
This property of a Hilbert-Schmidt operator L follows from the fact that L is 
the limit in the operator norm of finite rank operators. The space T-CS{Ei, E2) 
endowed with the Hilbert-Schmidt norm defines a Hilbert space. 

Some general references for this subsection are: [IE], [18], [30], [31] . 

2.2 Linear semi-groups and unbounded operators. 

Let L be a bounded linear operator on E. For every t G M, define: 

00 ^ 

$ (t) = e*^ = V -rt"i^", 
^ ^ ^ n! 

i=0 
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which converges in the operator norm. Then $ {t) is a bounded hnear oper- 
ator for every t, and we have the relation: 



$ (t + s) = $ (t) $ (s) Vs, t G M and $(0) = /, (7) 

where / is the identity operator on E, from which it follows that $ (t) and 
$ (s) commute and that $ (t) is invertible for every t. Relation (I7j) states 
that the map t i— >■ $ (t) is a group homomorphism. Note that it is continuous 
in the norm topology for operators: 

-/|| ^0 whent^O. (8) 

The solution of the Cauchy problem: 

^ = Lx{t), (9) 

x(0)=Xo (10) 

is given by x (t) = $(t)x(0). In other words, $ (t) is the flow associated 
with the ordinary differential equation Q. We can recover L from $ (t) by 
writing: 

Lx = \im-\^ (h) X - x] , x E E. (11) 

h^o h 

The norm continuity of the mapping t ^ $(t) is exceptional and has to be 
replaced by a more useful weaker property (cf. Definition 1, Sect. 1, Chap. 
IX of 



Definition 4 A family $ (t) ,t > 0, of bounded operators on E is called a 
one parameter semi-group if $ (0) = /, and for allt>0 and s > we have: 

$ (t + s) = $ (t) $ (s) = $ (s) $ (t) . (12) 

It is said to be strongly continuous or to be of class (Co) if, for every x E E, 
we have: 

lim$(t)x = x. (13) 

It is said to be a contraction semi-group if ||$(t)|| < 1 for all t > 0. 

Note that, since equality ( |T2i) is supposed to hold only for positive s and 
t, the operators $ (t) are no longer necessarily invertible, as in the case of 
a group. It can be proved easily that, if the semi-group $ (t) is strongly 
continuous, then lim^.^^ $ (s) x = $ (t) x and there are constants c and C 
such that 11$ (t)|| < Cexp (ct) . We also note that if [0, oo[ 9 t i-^ <l>(t) is a 
one parameter semi-group, so is the family of adjoint operators [0, oo[ 9 t t-^ 
$*(t), where we define $*(t) = ($(t))*. 
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Example 5 

In the situation of (i) (resp. of (ii)) of Example [H for given a, let $i(a) = A 
(resp. $2(0) = A). Then R 3 t 1— $i(t) is a strongly continuous contraction 
group. However [0, oo[ 9 t 1— >• $2(^) is only a strongly continuous contraction 
semi group, which can not be extended to a group. In fact, $2(^) is not 
invertible for t > 0. 

We now try to extend formula ( fTTl) . It turns out that when $ is no longer 
norm-continuous, but only strongly continuous, the right-hand side does not 
converge for every x, and if the limit exists, it does not depend continuously 
on X. The set of x for which the limit exists is obviously a linear subspace of 
E and on this subspace the limit is a linear function, let's say G of x. More 
formally, let 'D{G) be the subset of E of all elements x E E for which the 
strong limit 

Gx = lim - [$ (h) X - x] (14) 

exists. 

Theorem 6 Assume ^ is a strongly continuous semi-group. The set T>{G) 
is then a dense linear subspace of E and G given by [l4\ ) defines a linear 
map G : T>{G) E. This map is closed, i.e. if Xn is a sequence in T>{G) 
such that Xn X E E and Gxn ^ y E E then x G T>{G) and y = Gx. 
For every x G T>{G) and t > we have $ (t) x G V{G), 

G^ (t) X = ^ (t) Gx (15) 

and 

(t) X = G^ (t) X. (16) 

Proof. By definition T>{G) is the set of x where the limit in formula ffT^ 
exists (note that this is a strong limit, meaning that we should have norm- 
convergence), and Gx then is the value of that limit. Clearly G : 'D{G) — ^ E 
is a linear map. 

Given any x E E and t > 0, consider the integral: 

X{t)= [ <!>{s)xds. 
Jo 
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It is well-defined since the integrand is a continuous function from [0,t] into 
E. Using the semi-group property, we have: 



l[^{h)x{t)-x{t)] = l 



$ (/i) / $ (s) xds - $ (s) xds 
Jo Jo 

/ $ (s + /i) xds - $ (s) xds 
Jo Jo 

1 /"^ 

^{s + h)xds-- I ^ is) xds 



—y^{t)x-x. 

This proves that X (t) belongs to T>{G). Then so does jX{t), and when 
t — > 0, we have (t) — > X, so is dense in if, as announced. 

Now write: 



1 $ f/iV 



-X 



$(/.) 



h 



-$ (t) X. 



If X G V{G), the second term converges to $ (t) Gx and the third one to 
G$ (t) X. Formulas ([T^ and (fTUj) now follow, since these two terms must be 
equal. 

To prove the last condition, note that: 



VxgP(G'), <l>(t)x-x 



$ (s) Gxds. 



(17) 



Indeed, we have two functions of t, with values in E , which are zero for 
t = and which have the same derivative, namely $ (t) Gx, for every t > 0. 
So they must be equal. Now take a sequence x„ —>■ x, and assume that 
Gxn = Hn y E . Writing x = x„ in formula (fTTl) . we get: 



$ (t) X - X 



$ (s) yrfs. 



Dividing by t and letting t — > 0, we find that x G I^(G) and that y = Gx. 



Definition 7 In the situation of Theorem G is called the infinitesimal 
generator of the semi-group $. 

A linear map L : ViL) E2, where V (L) is a subspace of Ei, is called an 
operator from Ei to E2 with domain T> (L) . That two operators are equal, 
Li = L2, means that they have the same domain T>{Li) = T>{L2) and that 
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Lix = L2X for all x in the domain. The operator L is densely defined if 
T>{L) is dense in Ei. It is called a bounded operator if there exists a finite 
constant C > such that for all x G T^lL) one has ||i^a;|| < C||x|| and it 
is called an unbounded operator if such C does not exist. It is closed if its 
graph {{x, Lx) \ x E V (L)} is a closed subset of Ei x E2, which extends the 
definition in the preceding theorem. With these definitions, we can rephrase 
part of the preceding theorem by saying that every strongly continuous semi- 
group in E has a unique infinitesimal generator, which is a densely defined 
closed operator in E. The problem to determine if a given densely defined 
closed operator L in ii^ is the infinitesimal generator of a strongly contin- 
uous semi-group is more difficult and we refer the interested reader to the 
references mentioned in the end of this subsection. 

The definition of the adjoint of an operator can be extended to unbounded 
operators. Let L be a densely defined operator from Ei to £'2. We introduce 
the adjoint operator L* to L. The domain of V{L*) consists of all y E E2 for 
which the linear functional 

X (-^ (y, Lx) (18) 

is continuous on T>{L), endowed with the strong topology of Ei. For y G 
V{L*) we define L*y by 

{L*y, x) = (y, Lx) Vx G V{L). (19) 

This defines L*y uniquely, since T>{L) is dense in Ei. One proves that V{L*) 
is dense in E2 if L is also closed. 

An operator L in is called selfadjoint ii L* = L and skew-adjoint if 
L* = —L. We have the following clear-cut result (Stone's theorem): L is the 
infinitesimal generator of a group of unitary operators iff L is skew-adjoint. 

Example 8 

In the situation of Example [5l let Li and L2 be the infinitesimal generators 
of $1 and $2 respectively. Li is given by 

V{L,) = {/ G L\R) I f G L\R)}, 

and {Lif){x) = f'{x), where /' is the derivative of /. L2 is given by 

I?(L2) = {/gL2([0,oo[) If gL2([0,oo[)}, 

and {L2f){x) = f'{x). Since $1 is a group of unitary operators, we have that 
LI = —Li. $2 is not a semi-group of unitary operators, so L2 7^ —L2. A 
simple calculation shows that 

ViL*2) = {fe L\[0, oo[) I /(O) = and /' G L'i[0, oo[)} 
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and {Llf){x) = —f'{x). So here V{L*2) C T>{L2), with strict inclusion. One 
checks that $2 is a strongly continuous semi- group in L^([0, cxd[ ). It represents 
right translations of functions. Its infinitesimal generator is Lg. 

Some general references for this subsection are: [TB], [IE], [31], [S3] . 



2.3 Sobolev spaces 



For any integer n > 0, the Sobolev space if"(]R) is defined to be the set of 
functions / which are square-integrable together with all their derivatives of 
order up to n: 



k=l 



dx < 00. 



This is a linear space, and in fact a Hilbert space with norm given by: 

' dx^' 




k=l 



dx 



It is a standard fact that this norm of / can be expressed in terms of the 
Fourier transform / (appropriately normalized) of / by: 



2 



{i + yT fiy) 



dy. 



The advantage of that new definition is that it can be extended to non- 
integral and non-positive values. For any real number s, not necessarily an 
integer nor positive, we define the Sobolev space if'^(]R) to be the Hilbert 
space of functions associated with the following norm: 



{i + y'Y fiy) 



dy. 



(20) 



Clearly, H°(R) = L'^{R) and H'(R) C (M) for s > s' and in particular 
H'{R) C L2(M) c H^'Ir), for s > 0. H'{R) is, for general s G M, a space 
of (tempered) distributions. For example S^''\ the fc-th derivative of a delta 
Dirac distribution, is in H-^-^/^~'(R) for e > 0. 

In the case when s > 1/2, there are two classical results. 

Theorem 9 (Continuity of multiplication) Ifs > 1/2, iff and g belong 
to H'^(K), then fg belongs to iJ'^(R), and the map {f,g) fg from x H'^ 
to H'^ is continuous. 
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Denote by C^^(ffi) the space of n times continuously different iable real- valued 
functions which are bounded together with all their n first derivatives. Let 
C^q(M) the closed subspace of C^(M) of functions which converges to at 
±00 together with all their n first derivatives. These are Banach spaces for 
the norm: 

\\f\\cv=^^^ = max 0- 

b 0<k<n ^ ' ' 0<k<n " "^b 

Theorem 10 (Sobolev embedding) If s > n + 1/2 and if f E if*(]R), 
then there is a function g in C^q(M) which is equal to f almost everywhere. 
In addition, there is a constant Cg, depending only on s, such that: 

WgWc^ < c.II/IIh- 

From now on we shall no longer distinguish between / and g, that is, we 
shall always take the continuous representative of any function in iJ'^(]R). 
As a consequence of the Sobolev embedding theorem, if s > 1/2, then any 
function / in if'^(M) is continuous and bounded on the real line and converges 
to zero at ±cx), so that its value is defined everywhere. 

We define, for s e M, a continuous bilinear form on if~'^(M) x H''^(M.) by: 

<f,9>= r (Ky)) 9{y)dy, (21) 

where 1 is the complex conjugate of z. Schwarz inequality and ( l20l) give that 

\<f,9>\<\\f\\H-49\\H^^, (22) 

which indeed shows that the bilinear form in ( l2Ti) is continuous. We note 
that formally the bilinear form ( l2Ti) can be written 

/oo 
f{x)g{x)dx, 
-oo 

where, if s > 0, / is in a space of distributions H~^{M.) and (7 is in a space of 
"test functions" H'{R). 

Any continuous linear form g u{g) on H^{M.) is, due to (1201) . of the 
form u{g) =< / , (7 > for some / G if~^(M), with = so 

that henceforth we can identify the dual (i?^(M))' of H^(R) with H~'{R). In 
particular, if s > 1/2 then H'{R) C C^O^), so H~'{R) contains all bounded 
Radon measures. 

In the sequel, we will also be interested in functions defined only on the 
half-line [0,oo[. Let s > 0. We define the space i/'*([0,oo[) to be the set of 



14 



restrictions to [0, oo[ of functions in iJ^(M). This is clearly a linear space. To 
turn it into a Hilbert space, we have to use the following norm: 

ll/lk-Cp.oo)) = inf I gix) = fix) a.e. on [0, oo)}. (23) 

This is a Hilbert space norm on H^{[0, oo[), which is the natural restriction 
of the norm on if*(M). For instance, if / is a function in if*(]R) such that 
/ (x) = for X < 0, then its restriction /o to [0, oo[ belongs to H'^{[0, oo[), 
and we have: 

If s = n is an integer, the norm on H^{[0, oo[) turns out to be equivalent to 
the following one: 




To establish properties of translations in H^{[0, oo[ ), we need to know if there 
is a continuous linear embedding of iJ*([0, oo[) into iJ*(M), i.e. to know if 
the restriction operator has a continuous right-inverse. Fortunately, as we 
are in a Hilbert space setting, this problem is easy to solve. Let s > and 
let be the subset of functions in H^(M.) with support in ] — oo, 0], so that 
f e Ht_ if and only if / G H''(U.) and f{x) =0 for all x > 0. iff. is a closed 
subspace of Two functions /i, /2 G H^(R) have the same restriction 

to [0, oo[ iff /i — /2 G Hi. This means exactly that H^{[0, oo[) is a quotient 
space: -^'^([0, oo[) = H^(R) / H^. Introducing the notation © for the Hilbert 
space direct sum, we have the following result, which proof we omit since its 
trivial: 

Proposition 11 For s > we have: 
i) H'{R) = H'{[0,oo[)® Hi. 

a) Let M he the orthogonal complement of Ht in H''^{M.) w.r.t. the scalar 
product in H'^(M), let k be the canonical projection of H^(M.) on H^{[0,oo[) 
and let l be the canonical bisection o/ if^([0, cx)[ ) onto M. Then k is contin- 
uous, L is a Hilbert space isomorphism, kl is the identity map on -ff'^([0, oo[) 
and in is the orthogonal projection map in H'^iW) on M. 

We note that i is a continuous operator extending functions on [0,oo[ to 

functions on M and that ||/||jy.([o,oo[) = \\^f\\H-{R) ■ 

The dual space of -ff*([0, oo[) can easily be characterized in terms of 
distributions. For s > 0, H'{[0,oo[) = H'{R)/Hi, so 

{H-\[0,oo[)y={feH-\R)\ <f,g>=0\fgeH^_}. (24) 
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For s > 0, we define H ''([0, oo[ ) to be the closed subspace of all distributions 
in H~^{E?) with support in [0, oo[. It then follows that (if*([0, oo[ ))' can 
be identified with H~'{[Q,oo[). Since (i/^([0, oo[ ))" = i7^([0,oo[), it then 
follows that 

(i/"([0,oo[))' = i7-"([0,oo[) s gM. (25) 

If s G M, then the constant function taking the value 1 is not in -ff^([0, oo[). 
If s > 1/2, then even every function in i?^([0, oo[) converges to zero at oo. 
For this reason, we will need a larger class of distributions containing the 
constant functions. Let s G M and let / be a distribution with support in 
[0, oo[ such that it admits the decomposition f = g + a, where g G -^'^([0, oo[ ) 
and a G M. This decomposition of / is then unique and the set of all such 
distributions is naturally given the Hilbert space structure H^{[0,oo[) © M. 
The norm of / = + a is then given by 

11 r-\\2 II ||2 , 2 

11/11 — llfl'llH»([0,oo[) + " ■ 

This unique decomposition property leads us to the following 

Definition 12 For s G M, set E'^O, oo[) = i?^([0, oo[) © M with the corre- 
sponding Hilbert space norm. If f & i?'^([0, oo[) and if g & H''^{[0,oo[) and 
a G M are related by the unique decomposition f = g + a, then the norm of f 
is given by 

WfWls = \\g\fHs+a^. 

The dual (E^([0, cx)[ ))', of E^([0,oo[) is identified with (if^([0, oo[ ))' © M ^ 
£'"'^([0, oo[ ) by extending the bi-linear form, defined in (12T]) . to £'^''([0, oo[ ) x 

< F , G >= ab+ < f , g >, (26) 

where F = a + f e E-'{[0,oo[), G = b + g e £"([0,oo[), a, 6 G M, / G 
H-'{[0,oo[) and g G H'{[0,oo[). 

For all the Sobolev spaces H'^ we have introduced, and also for the spaces 
E''', there are two natural realizations of the dual space. Let us consider only 
the case of -E'^([0, cxd[ ), the other being similar. One possibility, the canonical 
one, is to identify (i?*([0, cxo[ ))' with E'^{[0,oo[) by the scalar product in 
i?*([0, oo[). This gives the Riesz representation in ([5]). Another possibility, is, 
as we have seen, to identify (i?'^([0, oo[))' with £'"''([0, oo[), by the bi-linear 
form defined in ( l26l) . There is a linear continuous map S : £''^([0, oo[) 
(£■^([0, oo[ ))' with continuous inverse, relating the two realizations. It is 
defined by: 

(/,^7)i^^ao,oo[) =<Sf,g >, \/f,g G £^([0,oo[). (27) 
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Now, different realizations of tlie dual space leads to different realizations of 
adjoint operators. Let A be a closed and densely defined operator from a 
Hilbert-space H to _E"*([0, oo[ ). We have already defined in ( fTSl) its adjoint 
operator A* from £'*([0,oo[) to H w.r.t. the duality defined by the scalar 
product. Let the dual H' of H be realized by Hi and the continuous bi-linear 
form < , >i: Hi x H ^ The adjoint A', w.r.t. the duality realized by 
< , >i and < , > is the operator from i?~'^([0, oo[ ) to Hi, defined by: 
The domain of V[A ) consists of all y G E~^{[0,oo[) for which the linear 
functional 

X <y, Ax > (28) 
is continuous on 'D{A). For y E 1^{A') we define A'y by 

<A'y,x >i=< y, Ax> \/x E V{A). (29) 

This defines A'y uniquely, since 'D(A) is dense in H. 

We now study translation semi-groups in the different spaces we have 
introduced. It follows directly from the definition fl2UI) of the norm in if*(M) 
and by dominated convergence that that left-translations defines a strongly 
continuous group of unitary operators £ in if^(M) for s G M (similarly as to 
the case of $1 in Example [5]): 

{jCtf){x) = f{x + t), V/ E H'{R) and t E R. (30) 

Since, for s > 0, the closed subspace Ht of H^{R) is invariant under the semi- 
group £i, t > 0, it defines a semi-group C in -ff*([0, oo[). Defining £ also on 
constants a G M by Ct^ = a we extend the semi-group C to £'^([0, oo[), 
s > : 

{Ctf)ix) = fix + t), V/GE^([0,oo[)andt>0. (31) 

Proposition 13 If s > 0, then C is a strongly continuous contraction semi- 
group on £'^([0, oo[). Its infinitesimal generator, denoted d, has domainV{d) = 
£■^+^([0, oo[). /// G £"'+^([0, 00 [) then df = /', where f is the derivative of 
/• 

Proof. We first observe that, in the canonical decomposition E^{[0,oo[) = 
H'^dO, oo[) ©M in Definition [T2l leaves the subspace H'^{[0, oc[) invariant 
and acts trivially on M. It is therefore sufficient to prove the statement with 
E"([0,oo[) replaced by iJ"([0,oo[). 

We use the notations of Proposition [TT] and let P = Lnhe the orthogonal 
projection on M. Since CtHi C i^i, for t > 0, it follows that PCt{I-P) = 0. 
The group composition law Ct'^u = ^^t+u, then gives for t, n > : 

{PjCtP){PCuP) = PCt+uP 
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So, [0, cxd[ 9 t (-^ PCtP is a semi-group of bounded operators on M. It is a 
strongly continuous contraction semi-group since this is the case for C and 
||P|| = 1. 

We have that Ct = K£ft, for t > 0. Using that Ct = uPCtPi it eas- 
ily follows from the semi-group property of PCtP that £ is a semi-group 
on if'^([0, oo[ ). It is a strongly continuous contraction semi-group, since 
this is the case for PCP and since ||k|| = = 1. Let d be the in- 
finitesimal generator of C. By the definition of C it follows that T>{d) = 
{/ e H%[Q,oo[) I f e H'{[0,oo[)} and df = f for / G V{d). But 
H'+\[0,oo[) = {f e H'{[0,oo[) I /' G if"([0,oo[)}, which proves the propo- 
sition. ■ 

Example 14 

Let C'^ : E-^([0,oo[) ^ E-"([0,oo[) be the adjoint of Ct, t > in Propo- 
sition [131 w.r.t. duality defined by the bilinear form < , > . £' is then a 
semi-group of right-translations on the space of distributions £'"''([0, oo[ ). 
Loosely speaking {C[f){x) = f{x — t). Let s > 1. hen the generator d' 
has domain i?~'^"'"^([0, oo[ ) and —d' is the derivative of distributions, so 
{& f){x) = —df{x)/dx if / is a differentiable function. One is easily con- 
vinced that the expressions for CI and d* are more complicated. 

Some general references for this subsection are: PP, [S], [13] • 

2.4 Infinite-dimensional Brownian motion 

In this sub-section we consider a separable Hilbert space E and an index- 
set I with the cardinality equal to the dimension of E. The space E can be 
infinite-dimensional or finite-dimensional. There is given a family W^, ^ G I of 
standard independent Brownian motions on a complete filtered probability 
space (fi, P, JF, ^). The filtration A = {J^t}o<t<T, is generated by the W\ 
and J-' = Tt- 

Definition 15 ^4 standard cylindrical Brownian motion Wt, < t < T, 
on E is a sequence CiW^, z G I o/ E-valued processes, where the Cj are 
the elements of an orthonormal basis of E and the W^, z G I, are indepen- 
dent real-valued standard Brownian motions on a filtered probability space 

in,p,T, A). 

From now on, given a standard cylindrical Brownian motion W, we shall 
write informally Wt = Xliei ^t^i- ^ is finite, we have: 

\\Wtf = J2\\^t\\'^ < °° 
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and Wt is a stochastic process with values in E. If I is infinite, then for every 
t the right-hand side is the sum of infinitely many i.i.d. positive random 
variables, which does not converge in any reasonable way. In that case, the 
formula Wt = X^iei ^t^i cannot be understood as an equality in E, and must 
be given another meaning. 

Proposition 16 // Wt = ^jgj W^ei is a standard cylindrical Brownian mo- 
tion, then, for every f E E with \\f\\ = 1, the real-valued stochastic process 
W/ defined by 

W/ = J2ie^J)Wi (32) 
is a standard Brownian motion on the real line. 

Proof. If I is finite, the result is obvious. Let us then consider the case 
when I = N. We first have to check if the right-hand side is well-defined. By 
Doob's inequality for martingales: 



E 



sup 

0<t<T 



n+p 



< AE 



n+p 



J2{e.J)W^ 

i=n 
n+p 

<4Tj](e„/f ^0. 



This implies that the right-hand side of (1321) converges in probability to a 
continuous process. Since each finite sum is Gaussian, so is the limit, and 
the result follows. ■ 



So, in the case when I is infinite, the r.h.s. of Wt = X^jgj W^Ci makes no 
sense in E, but every projection does. Equation flH^ can be rewritten as: 



V/eE, {WtJ) = Y,{e.J)Wi 



We will now show that the stochastic integrals with respect to cylindri- 
cal Brownian motion make sense, provided the integrand satisfies a strong 
integrability condition. Consider the space T-CS{E,F) of all Hilbert-Schmidt 
operators from E into a Hilbert space F. Let the space {T-[S{E, F)) consist 
of all progressively measurable processes A with values in the Hilbert space 
nS{E,F), such that: 



E 



T 



\A 



tWns 



dt 



< oo. 



19 



Recall that we have, according to Definition [SI 



I A 



t\\HS 



n=0 



t^nW , 



where (cr, 



'nm 



is any orthonormal basis of E. 



Theorem 17 The stochastic integral: 

cT 



AtdWt 



is well-defined for every process A E {T-[S{E, F)) . It is a continuous mar- 
tingale with values in E, and we have the usual isometry: 







[ AtdWt 




Jo 


L2 ^0 



E[\\At\\i,s\ dt. 



In other words, the random variable AtdWt has mean and its variance is 

X]^o ^ [II Ate„||^] dt, the sum of the variances of the independent sources 
of Gaussian noise. 

As usual, by localization the stochastic integral can be extended to a 
wider class of processes. Denote by Cf^^ i'HS) the set of all progressively 
measurable processes with values in TiS, such that: 



P 



I^Whs dt < oo 



1. 



Then the stochastic integral defines a continuous local martingale. 
Some general references for this subsection are: [7j, [H], [23] 



[26]. 



3 The dynamics of bond prices 
3.1 The non- parametric framework 

From now on, and for the rest of the paper, we are given a finite time interval 
of possible trading times T = [0,T] and we are given a family W^, i G I of 
standard independent Brownian motions on a complete filtered probability 
space {Q, P, JF, A), the filtration A = {^t}o<t<T^ is generated by the W\ and 
JF = J^f. The family I itself can be finite or infinite, in which case we take 
I = N. Let -^^(I), be the Hilbert space of all real sequences x = {xi)i^i, such 
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that ||x||f2(-j) = C^i(zi{o-nYY^'^ < oo. So, when I has a finite number m of 
elements, then £^{T) = M"^. Often we write just £^ for 

Heath, Jarrow and Morton (henceforth HJM) were the first to study the 
term structure of interest rates in a non-parametric framework. Their basic 
idea (see [I2]) consists of writing one equation for the price of every zero- 
coupon at time t. Denoting by BtiT) the price at time t of a zero-coupon 
bond maturing at time T > t, the HJM equation has the following form: 

Bt{T)=Bo{T)+ [ Bs{T)as{T)ds+ [ Y,Bs{T)vl{T)dWl 0<t<T 
Jo Jo .^j 

(33) 

There are infinitely many such equations, one for each maturity T >t. 

The trend at (T) and the volatilities vl(T) are supposed to be progressively 
measurable processes, which means, for instance, that they could be functions 
of all the Bs{S), S > s and s < t. In due course, we will make further 
assumptions so as to ensure that equations such as (l33l) make mathematical 
sense. 

Let us discount all prices to t = 0, by the spot interest rate rj, which in 
terms of the zero-coupon bond price is given by 



OBtiT) 

rt = 



dT 

The discounted prices of zero-coupons are now: 

ft 



(34) 

T=t 



Bt{T) = Bt{T)exp{- [ rjs) (35) 

Jo 

and the equations (IHHjl become: 



Bt{T) = Bo{T)+ [ B,{T){as{T)-rs)ds+ [ Y,B,{T)vi{T)dWl 0<t<T 
Jo Jo .gj 

(36) 

and, again, there is one such equation for every maturity T > t. Note the 
boundary condition Bt(T) = 1, and hence, from ( l35l) : 

Bt{t) = exp(- / rjs). (37) 



3.2 The bond dynamics in the moving frame 

For every x > 0, we denote by pt{x) the price and by pt{x) the discounted 
price at time t of a zero-coupon maturing at time t + x. The stochastic 
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processes Bt{T) and pt{x) are related by: 

Pt{x) = Btit + x). 

In other words, as explained in the introduction, instead of dating events by 
their distance from a fixed origin, defined to be t = 0, we are dating them by 
their distance from today: we are using a time frame which moves with the 
observer. The equation for pt in the moving frame, is easily obtained from 
(l36l) . For every x > 0, we have: 



Pt{x) =po{t + x) + i ps(t - s + x)ms(t - s + x)ds 
Jo 

+ / ^p,(t-s + x)a^(t-s + x)c?iy;. 



(3^ 



where 

mt{x) = a{t,t + x) — rt and al{x) = v\t,t + x), (39) 

for all < t < T and x > 0. Here, again, the trends t i— > mt{x) and the 
volatilities t crl{x) are progressively measurable processes. 

Instead of looking at ( l38i) as an infinite family of coupled equations, one 
for each x > 0, we shall interpret it as a single equation describing the 
dynamics of an infinite-dimensional object, the curve x i— ^ Pt{x), which will 
be seen as a vector pt in the Hilbert space £''^([0, cx)[ ), for some fixed s > 1/2, 
chosen so that the functions rrit and al belong to E'^{[0,oo[). 

Let C be the semi-group left translations on E^{[0,oo[) (see formula 
(l3T]l and Proposition [T3ll . From now on we shall just wright E'^ instead 
of £^^([0, oo[ ), when there is no risk of confusion. The equations in ( l38l) can 
be rewritten as one equation in E^: 

Pt = CtPo+ [ {Ct.s{Psms))ds + [ J2(^t-s{Ps(rl))dW:. (40) 
Jo Jo .gj 

Theorem 18 Let s > 1/2. Assume that po G E^ and assume that rrit and 
the (Tf, i G I, are progressively measurable processes in E'^ satisfying: 

rf 

I {\\mt\\E^ + ^\\(yl\Ws)dt < (x> a.s. (41) 
Then equation defines a unique process p in E" satisfying: 

/ + \\Ptmt\\Es + V \\pt(JlfEs)dt < oo a.s. (42) 

"^0 m 
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The process p has continuous trajectories in E'*, 




Pt = exp{ / Ct-s {{ms-- y(aiy)ds + V aidW^ } C^Po- (43) 



and if Po G then the process p takes its values in H'^. If po G satisfies 
> (resp. Po > 0), i.e. po{x) > (resp. po{x) > 0) for all a; > 0, then so 
does pt- 

For a proof of this theorem see Lemma A.l of [10], which is reproduced in 
the appendix of this article (Lemma HH]). Note that equation (l40l) imphes 
that Po is the value of pt for t = 0. 

A word here about the choice of function spaces. Assuming that pt belongs 
to H'^ for some s > 1/2 is minimal: it is basically saying that the zero-coupon 
prices depend continuously on time to maturity and go to zero at infinity. 
This, however, is too strong a requirement for rrit and the a^: we cannot 
expect the trend and the volatilities to go to zero when the time to maturity 
increases to infinity. This is why we are assuming that and the belong 
to E'^. To simplify the mathematical formalism and also to include interest 
rate models, with vanishing long term rates, we have permitted that pt & E''. 
Now according to Theorem [181 Pt is in-fact in if* if po G H'^ . 

Condition (HTI) implies that Xltei ll^tlls^ finite for almost every (t,uj) G 
T X fi. This means, when I = N, that the operator cr^ from £^(I) to E'^ defined 
by: 

atCi = al z G I, (44) 

where Cj are the elements of the standard basis of ^^(I), is Hilbert-Schmidt 
a.e. {t, Lj). We have 

11 ||2 _ II i||2 

\\"t\\HS{£^,E'') ~ 2^ Wt\\E^- 

We shall refer to a as the volatility operator process. It takes its values in 
T-CS{i'^, E'^) and when we say that it is progressively measurable, it is meant 
that all the cr* are progressively measurable. 

We can now, using the stochastic integral introduced in Theorem [TTl 
rewrite equation (l40i) on a more compact form in E'^ , where s > 1 /2 : 



Pt = CtPo+ / Ct-s{Psms)ds + / Ct-s{PsO-s)dWs. (45) 
Jo Jo 

This makes sens in Indeed, the only difference with equation fl4Up is the 
last term on the r.h.s. When condition ( HTl) is satisfied then the volatility 
operator au, defined by (jSj), from £^ to E'^, is Hilbert-Schmidt a.e. {u,uj). 
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Since pointwise multiplication of functions in is a continuous operation 
for s > 1/2 it follows that the linear operator x ^ PuCr^x, from i"^ to -E*, is 
Hilbert-Schmidt a.e. (m, uj). L„ is bonded for every f > 0, so the integrand is a 
progressively measurable l-iSilP ^ £''^)-valued process satisfying the conditions 
of Theorem [T71 

A process p with values in E'^ satisfying (1451) (or equivalently (1401) ) and 
(I42p will be called a mild solution of the bonds dynamics. 

Note that we are not worrying about the boundary condition ( !37j) at this 
time, because it does not make mathematical sense: how do we define rj ? 
This will be taken care of in the next section. 

3.3 Smoothness of the zero-coupon curve. 

Another way to proceed is to write ( l38l) in differentiated form. For fixed 
X > 0, a formal calculation using Ito's lemma and which can be rigorously 
justified gives: 

dpt{x)-pt {x) mt (x) dt - y^/pt{x)al{x)dWl 



In the expression on r.h.s. we can replace d/dt by d/dx, since po and the 
integrands on the r.h.s. are functions oi t + x. Derivation w.r.t. x under the 
integral then gives: 



The l.h.s. is equal to {{d / dx)pt{x)) dt, according to (l38l) . so 

d 

dpt{x) - pt (x) mt (x) dt - ^pt{x)al{x)dWl = {^Pt{x)) dt, (46) 
for all X > and t G T. 
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Introducing the infinitesimal generator d of tlie semi-group C (see Propo- 
sition [T3]), tliis can be understood as an equation in E'^ : 

dpt = (dpt + Ptmt)dt + J2 Pt^yldWi (47) 

or equivalently: 



Pt=Po+ / {dps + Psms)ds + / ^PsaldW'^. 
Jo Jo 



(48) 



Equation (HUI) is tlie integrated version of (HHj) . w.r.t. the semi-group C. 
The connection between formulas (HHl) and (HOl) is similar to the variations 
of constants formula for ODE's in finite dimension. 

We now have to give some mathematical meaning to equation ( HHl) . This 
will require beefing up the existence conditions given in Theorem [THJ The 
following corollary follows from applying Theorem [TH] with s + 1 instead of 
s : 

Corollary 19 Let s > 1/2. Assume that p^ G 'D{d) = E^~^^ and assume that 
rrit and the cr^, i G I are progressively measurable processes with values in 
E''^'^^ satisfying 

/ (||"2t||ijs+i + ^ \\al\\%+i)dt < oo a.s. (49) 

Then the mild solution p, in Theorem\TE, of the bonds dynamics satisfies the 
following condition: 

pt G E'+^ and / (|bi|U»+i + \\ptmt\\E^ + IbtC^tlllO dt < ^ a.s. (50) 

Equation (^5]) holds for every t. In addition p has continuous paths in E'^~^^ 
andpt G H^~^^ if Pq G 

By definition a solution of equation fj^Uj) is called a strong solution of the 
equation fHHjl . when condition fl3Ul) is satisfied. Here we shall say that p is a 
strong solution of the bonds dynamics. 

As a consequence, in the situation of Corollary [191 the term structure x i— > 
Pt{x) is for every t, and interest rates are well defined. The instantaneous 
forward rate Rt{x) contracted at t G T for time to maturity x and the spot 
rate at time t, for instance, are defined by: 

dlogptjx) idpt){x) (dpt) (0) 

ox pt{x) pt{0) 
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By Corollary (THl p is a strong solution and the maps t ^-^ pt and t dpt 
are continuous from T into E^, and hence into C'^([0, oo[) endowed with the 
topology of uniform convergence. So Ps{0) and (dps) (0) converge to pt{0) 
and (dpt) (0) , when s ^ t. In other words, is a continuous function of t, 
when pt{0) > for all t e M. 

We are now able to make sense of the boundary condition (1H7|) . which we 
rewrite in terms of p : 

Pt{0) = eMf^-^^ds), (52) 

Jo Ps\^) 

for every t G T. 

Proposition 20 Let s > 1/2. Assume that rrit and the al are progressively 
measurable processes with values in E'^'^^ satisfying [J^ ) ^'^^ 

mtiO) = 0, ^^'(O) = \/iEl (53) 

and assume that Pq satisfies 

Po G E'+\ po(0) = 1, Poix) > Vx > 0. (54) 

Then the solution of the bond dynamics, given by Corollary [W[ satisfies the 
boundary condition [U^) . 

Proof. Since and the al take values in i?*'^^, they are continuous func- 
tion on [0,oo[, and condition 0531) makes sense. As po > it follows from 
Proposition [IH] that pt > 0. We have shown that, if pt is a strong and strictly 
positive solution of the bond dynamics, then vt given by ( 15T1) is a continuous 
function of t. Writing conditions (l53il into equation (HHl) . we get: 

Pt(0)=Po(0)+ f {{dp,m+pMmsmds+ [ Vp,(0)a:(0)diy; 
Jo Jo .gj 

= 1+ [ idps)iO)ds = 1 - I r,p,iO)ds. 
Jo Jo 

In other words, (p{t) = pt{0) must satisfy the differential equation (f'{t) = 
—rtif(t), with the initial condition (p{0) = 1. The result follows. ■ 

When we get to optimizing portfolios, we will need estimates on the 
solutions of the bond dynamics. They are provided by the following result: 
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Theorem 21 Let q{t) = pt/CtPo and q{t) = pt/CtPo. If Po, cr and m in 
Proposition [10 also satisfy the following additional conditions: 



E{{ [ \Wt\\ns(i^Es+i)dty + exp(a [ ||cTj^5(^2,Es)(it)) < cx), Va G [1, oo[ 
Jo Jo 

(55) 

and 

E{{ i ||m( ||£;s+i(it)" + exp(a / ||m( H^sc/t)) < oo, Va G [1, oo[ , (56) 
Jo Jo 

then the solution p in Proposition \^ has the following property: 

p, p, q, q, 1/q, 1/q G L^{n, P, L^{T, E'^')),Wu G [1, oo[ . (57) 
Proof. We use the notation 

St{L) = exp( r A-.((m, - ^ J](a^)2)rfs + aJW,)), (58) 

>'0 -^IT 



for 

-t 



Lt= f {rusds + aJW^), if < t < T. (59) 
Jo 

Conditions {i) — {iv) of Lemma HH] are satisfied for p. Estimate 01451) of Lemma 
HHlthen shows that pe L"(fi, P, L°°(T, E"+^)) Vm G [1,oo[. By the exphcit 
expression ( H3l) . q = S{L), so it follows from Lemma 149) that the conclusion 
holds true also for q. 

Let N, = /J((-m, + - E^el'^ldW:). Then 1/q = £~(iV). 

According to conditions (1^^ . (1561) . the conditions (i) — (iv) of Lemma HH] 
(with instead of L) are satisfied. We now apply estimate 01451) to 1/q, 
which proves that 1/q G L"(l], P, L°°(T, ^'+^)), for all u > 1. 

To prove the cases of g", a = 1 or a = —1, we note that q(t) = q(t)pt{0). 
Using that the case of is already proved and Holders inequality, it is 
enough to prove that g G L"(f2, P, L°°(T, M)), where g{t) = {pt{0))~". Since 
pt{0) = (Apo)(0)(g(t))(0) = po{t){q{t)){0), it follows that 

0<g{t) = {po{t))--{{q{t)m)-''. 

By Sobolev embedding, pq is a continuous real valued function on [0, oo[ and 
it is also strictly positive, so the function t i-^ (po(^))~" is bounded on T. 
Once more by Sobolev embedding, ((g(^))(0))~" < C|| (g(t))~"||£;s. The result 
now follows, since we have already proved the case of The case of p is so 
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similar to the previous cases that we omit it. ■ 

Under the hypotheses of Proposition [20l pt{0) satisfies ( l52|) . so it is the 
discount factor ( 1371) . It has nice properties, as follows from the second part 
of the proof of Theorem [2T] 

Corollary 22 Under the hypotheses of Theorem [21[ z/a G M, then the dis- 
count factor pt{0) satisfies 

E(sup(pt(0))°) < oo. 

iGT 

Remark 23 It follows from Theorem that for all t & T, pt and Pq have 
similar asymptotic behavior. In fact for some r.v. A > 0, A-^po{t + x) < 
Pt{x) < Apo{t + x) , for allt and x > 0, where A is independent of x and 
t and Ae L^'iVL, P) for all u>l. 

In a different context, Hilbert spaces of forward rate curves were consid- 
ered in [1] and [H]. The space E"^, with s > 1/2 sufficiently small, contains 
the image of these spaces, under the nonlinear map of forward rates to zero- 
coupons prices. Or more precisely, it contains the image of subsets of forward 
rate curves / with positive long term interest rate, i.e. f{x) > for all x 
sufficiently big. 

4 Portfolio theory 

In this section s > 1/2, = E"[0,oo)[) and T = [0, f], where f is the 
time horizon of the model. We also write E for E'"^ = -E''^[0, oo)[) and E' for 
E-^[0,oo)[). 

4.1 Basic definitions. 

We recall that, by the bilinear form < , >, the space E~^ is identified with 
the dual of E''', that is, the space of continuous linear functionals on E'''. 
It is important to note that, since s > 1/2, the space E'^ is contained in 
C°([0, cx))[), the space of bounded continuous functions on [0,oo[, so that 
E~^ contains the dual of C°([0, oo[), which is the space of bounded Radon 
measure on [0,cx)). In particular, all Dirac masses S^, for x > 0, belong to 
E-\ 
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Definition 24 A portfolio is progressively measurable process on the time 
interval T, with values in E~'^. If 9 is a portfolio, then its discounted value 
at time t is 

Vtie)=<9t,pt> . (60) 
The basic example is a portfolio of one zero-coupon: 
Example 25 

Consider a portfolio containing exactly one zero-coupon bond with maturity 
date T, i.e. time of maturity T : 

1) Let T >T and let T be fixed. The portfolio 6 is then defined by 

9t = 5t-u Vt < f . (61) 

Since T > T, we have indeed that the support of the distribution 6t is 
contained in [0, oo[, so 6t G -E'"'^. With this definition, the value of the zero- 
coupon is: 

< ST-t,Pt > = Pt{T -t) 
which is precisely what we had in mind. 

2) Let T < T and let T fixed. In this case we note that the process in ( 1611) 
does not continue after time T: the zero-coupon is converted into cash. So 
the buy-and-hold strategy is not possible for zero-coupon bonds, unless the 
horizon T is less than the maturity T. 

3) Let T = t + X and x > a fixed time to maturity. Then the portfolio is 
defined by 

Ot = S^, for t < f. (62) 

We note that the higher we choose s, the more portfolios can be incor- 
porated into the model. For instance, if s > 3/2, all curves in are C^, so 
that the derivative 6'^ of the Dirac mass belongs to E~^. The value of 
is: 

< 5'j._t,pt > = p[{T -t) = -Rt{T - t)pt{T - t), (63) 

where p'^ix) = dpt{x)/dx and where Rt{x), defined in (EI]), is the instanta- 
neous forward rate with time to maturity x, contracted at time t. This also 
implies that the higher we choose s, the more interest rates derivatives can 
be incorporated into the model. If s > 1/2, then we can contract directly on 
the values of zero-coupon bond prices, and if s > 3/2, then we can contract 
directly on the values of interest rates. 

We next introduce the notion of self-financing portfolio. We state a defi- 
nition such that it will makes sense for mild solutions of the bonds dynamics: 
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Definition 26 A portfolio is called self-financing if, for every t G T 



Vtie) = Voie) + [ <es, p.m. ds + J^Ps^dW: > . (64) 

Given a strong solution p of tlie bonds dynamics, we have for a self-financing 
portfolio: 

dVt{9) = <9t,dpt- dpt dt > . (65) 

Note that this is not the standard definition: this is because we are in the 
moving frame. Changes in portfolio value are due to two causes: changes in 
prices, as in the fixed frame, and also to changes in time to maturity. 

For the right-hand side of 0641) to make mathematical sense and to intro- 
duce later arbitrage free markets, we need a further definition. 

Definition 27 A portfolio 9 is an admissible portfolio if ||^||p < oo, where 



l = E 



[[ \<9t, ptmt > \dtY + [ V(< 9t , ptal >fdt 
Jo Jo .^j 



P is the linear space of all admissible portfolios and Pgf the subspace of self- 
financing portfolios. 



The discounted gains process G, defined by 



G{t, 9) = [ {<9s, Psm, > ds+ < 9^ , PsCr.dW^ >), 
Jo 



(66) 



is well-defined for admissible portfolios: 



Proposition 28 Assume that po, m and a are as in Proposition EO- If 
6* G P, then G{-,9) is continuous a.s. and E{sup^izj{G{t,9))'^) < oo. 



Proof. Let 6^ G P and introduce X = sup^^j \G{t, 9)\, Y{t) = < 9s , Psfng > 
ds and Z{t) = /J < 9^ , P.cr.diy, > . Then G{t,9) = Y{t) + Z{t), according 
to formula ( 1661) . Let p be given by Proposition [201 of which the hypotheses 
are satisfied. 

We shall give estimates for Y and Z. By the definition of P : 

E{{snp{Y{t)f) < E{{ r \<9s, Psfn, > \dsf) < \\9\\l. (67) 

tGJ Jo 
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By isometry we obtain 

Jo 



E{ fy,{< Os,PsK >)'ds)<\\e\\l. 



Doob's inequality and inequality (lUHl) give E{sup^^j Z{ty) < 4||^||p. In- 
equality (j671) then gives E^X"^) < 10||6'||p, which proves the proposition. ■ 



Example 29 

1) The portfolio in 1) of Example 1291 is self-financing and the portfolios in 2) 
and 3) of Example [29] are not self-financing. 

2) The interest rate portfolio in formula (l63l) is self-financing. 



4.2 Rollovers 

Definition 30 Let S > 0. A S-rollover is a self-financing portfolio 6 of 
a number of zero-coupon bonds with constant time to maturity S and with 
initial price Vq{6) = Po{S). 

It follows directly from the definition that a S'-rollover have the same initial 
price as a zero-coupon with maturity date 5*. It also follows that, if Xt is the 
number of zero-coupon bonds in the portfolio at t, then we must have: 

Ot = Xt6s, 

where the real-valued process x makes the portfolio self-financing. 
Proposition 31 If 6t is a S- rollover, then: 

Xt = exp( [ Rs{S) ds). (69) 
Jo 

Proof. The portfolio 6t only contains zero-coupons with time to maturity S, 
so that Vt{9) = XtPt{S). Assuming the process x to be of bounded variation 
it follows that: 

dVt{9)=ptiS)dxt + XtdptiS). 
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Substituting the expression for dpt{S) this becomes: 



dVtie) = Pt{S)^dt + xtd^pt{S)dt + xtPt{S){mt{S)dt + ^ (yl{S)dWl) 

= (PtiS)^ + xAptiS) + xMS)mt{S))dt + XtPtiS) J2 ^liS)dWi. 

According to flM|l the portfoho is then self-financing if and only if: 

dxt 



This means that: 



Pt{S)^+Xt{dp,){S)=0. 



1 dxt _ 1 dpt{S) _ ^^^^^ 



Xt dt Pt{S) OS 

and the formula ( l69i) follows by integration. This proves the proposition since 
X then is of bounded variation. ■ 

In particular, if S" = 0, then we get the usual bank account with spot rate r^. 

Henceforth, we will denote by qt{S) the value (discounted to t = 0) at 
time t of a S-rollover. In the preceding notation, qt{S) = Vt{6). 

Introducing the price curve of the roll-over at time t, qt : [0, oo[^ M, we 
find that the price dynamics of roll-overs is given by: 

qt=Po+ [ qsrrisds + I qs^ (^IdWl, (70) 
Jo Jo 

Note that, compared to the same formula for bond prices, the term in d has 
disappeared from the right-hand side. 

A S'-roUover is a bank account which needs advance notice to be cashed: 
if notice is given at time t, the rollover will then pay Xt units of account 
at time t + S. In other words, at time t, when notice is given, the rollover 
is exchanged for qt{S)/pt{S) = Xt units of a unit zero-coupon with time of 
maturity t + S. 

As we noted earlier, zero-coupons do not in general allow buy-and-hold 
strategies. However rollovers do: a constant portfolio of rollovers is always 
self-financing. A general bond portfolio 9t can be expressed in terms of a 
portfolio of rollovers rjt and vice versa. 
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4.3 Absence of arbitrage opportunities. 

Let p be a mild solution of the price dynamics. Suppose that 6t is a self- 
financing portfolio such that, for almost every (t, cu) G T x f2, we have: 

Vz G I, < Ot {u), Pt{uj)<jl{uj) >= 0. (71) 

(We note that Ptiuj) G -E^ is a function of time to maturity, x i— > pt{uj,x), 
and similarly for 9t etc.) Then f lM|l gives dVt{9) =< 6t , rritPt > dt, so that 
6t is risk-free. Since the spot rate is zero (after discounting values to t = 0), 
in an arbitrage free market it must follow that for almost every {t,uj): 

< Otiu) , mt{u)pt{uj) >= 0. (72) 

Comparing flTTl) and fl72|) . we find that pt{u)mt{uj) must belong to the closure 
of the linear span of {pt{u})al{uj) | z G I}. In fact this follows rigorously using 
Lemma [Ml proved independently of this subsection. There are now two 
cases: 

• I is finite. Then the linear span is finite-dimensional, and it coincides 
with its closure. So there are numbers '-fl{uj),i G I such that 

Pt{uj)mt{uj) = pt{uj) '^'^1{uj)(tI{uj) (finite sum). 
Since Ptifjj) > for almost every (t,^), this leads to: 

and since the processes m and a* are progressively measurable, so can 
one choose the processes 7*. Note that the preceding equation holds in 
E^, and that it translates into a family of equations in [0, oo[: 

mf(a;, x) = 7j(ci;)(tJ(co', x) Vx > 

or, as usual, omitting to mention the uo variable: 
mt(x) = ^7jVJ(x) Vx > 0. 

The 7^ are the components of a market price of risk, and they do not 
depend on the time to maturity x. Using the volatility operator process 
a the last equality reads 

mt = atit Vt G T (73) 

and any 7, progressively measurable with values in £^(I), satisfying this 
equation is called a market price of risk process. 
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• I = N. Then the hnear span is not closed in general; in fact, it is closed 
if and only if it is finite-dimensional. In that case, we shall impose a 
stronger condition. To prove that the market is arbitrage-free, we shall 
use that mt{uj) is in the range of the volatility operator at{uj) which is 
a subset of the above closed linear span. So, once more we impose that 
the condition (173|1 should be satisfied, but for 7 with values in If 
the range of 17^(0;) is infinite dimensional, then this condition is indeed 
stronger, since at{w) is a.e. a compact operator. 

In both cases, we also need that 7 satisfy some integrability condition in 
{u!,t). This leads us to the following 

Definition 32 We shall say that the market is strongly arbitrage-free if there 
exists a progressively measurable process 7 with values in such that 

mt = aau VteT (74) 

and 



E 



T 

2 



exp(a / ||7t||£2(it) 







< 00, Va > 0. (75) 



If the market is strongly arbitrage-free then, by the Girsanov theorem, a 
martingale measure is given by dQ = ^fdP, with: 

it = exp I -\ f hsWlds - > l^ldWl 1 . (76) 



2 




The W\ i e I, where 



Wi = W^+[\ids, (77) 

are independent Wiener process with respect to Q. The expected value of a 
random variable X with respect to Q is given by: 

Eq[X] = E[^tX]. 

Under a martingale measure, the discounted zero-coupon price process p 
satisfies the equation 

Pt = CtPo+ [ Ct^siPs(Ts)dW, (78) 
Jo 

and also the equation 

Pt=Po+ / dpsds + / psas dW^. (79) 
Jo Jo 
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The discounted roll-over price process qt is given by: 



qt=Po+ / qs(TsdWs. (80) 



Lemma 33 A portfolio 9 is self-financing if and only if: 

Vt{e) = v^{e) + / V < , vt<yl > dw:. (8i) 

We note that the integrand is in fact the adjoint operator of the operator 
bt{uj) = pt{uj)at{uj) from ^^(i) to E' {[0 , oo[) : 

{btiujyetY =<et,ptai>, Vzgl (82) 

To see this, with xl{uj) =< 6t , ptcrl >, rewrite it as follows: 
for all {t,uj) and all z G ^^(I) 

{z, Xt{uj))^2 = ^z' <et (lo) , pt (uj) o\ {yj) > 
=< Ot (cu) , pt {oj) ^ ai {uj) z' > 

=< Ot (uj) , bt (uj) z >=< bt (uj)' Otiuj), z > . 

If the market is strongly arbitrage- free and if condition ( 1551) of Theorem [2T] 
is satisfied, then also condition fl3Ul) is satisfied and the Theorem HT] applies. 

5 Hedging of interest derivatives 

From now on, it will be a standing assumption that po satisfies condition 
(|5^ . that a satisfy conditions fl53l) and fl55|l and that the market is strongly 



arbitrage-free according to Definition | 

Before we solve the optimal portfolio problem, we shall study the problem 
of hedging a European interest rates derivative with payoff X at maturity 
T. X is said to be an attainable contingent claim or derivative if Vf{6) = X 
for some admissible self-financing portfolio 6. Here we are only interested 
in payoffs, relevant for the optimal portfolio problem considered in these 
notes, i.e. X G L^(r2,jF, P) for every p > 1 (see Lemma HTj) . We first 
introduce the hedging equation, the Malliavin derivative and the Clark-Ocone 
representation formula, which then permits the reader, if he wish, to proceed 
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directly to the study of the optimization problem in the case of deterministic 
a and 7 in §6.2.11 

Assume that X G L'^{Q,T,Q), where Q is one equivalent martingale 
measure given by ( 176|) . Then, by the martingale representation theorem, X 
can be written stochastic integral: 

X = Eq[X]+ fY.'^.dWl (83) 

with: 

Eq[ / WxtWldt] < 00. (84) 
Jo 

Comparing with equations ( IHTj) and ( l82l) for a self-financing portfolio, we 
obtain the hedging equation 

ht{uj)'et{uj) = xt{uj), a.e. {t,u}), (85) 

where the operator bt{uj) = pt{uj)cr^{uj) from £^(I) to £'"([0, oo[ ) was intro- 
duced in (!82l) . Equivalently: for almost every (t, cu), 

xj(u;) = < (a;) , pt (u) a\ (a;) >, Vi G I . 

We next introduce the Malliavin derivative (c.f. [26]), D^X^ with respect 
to VT, at time t G T of certain T = Tf measurable real random variables X 
by: 

Dl) DtX = 0, if X is a constant, 

D2) DtX = hu like L2(T, ^{1)) and X = KdWl, 
D3) Dt{XY) = XDtY + YDtX. 

The algebra of such random variables is dense in L'^{Q,J^,Q), which can be 
used to extend the definition to larger sets. DtX takes its values in £^(I). 
The partial derivative, with respect to W^, Di tX, is the i-th component of 
DtX. 

We will use the following expression for the Malliavin derivative of an Ito 
stochastic integral: 

Dt r J2 <dW: = xt+ r J2(Dtxl)dW:, (86) 
when almost all the xl are Malliavin differentiable and sufficiently integrable. 
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In the case when X is Malhavin different iable, the Clark-Ocone represen- 
tation formula states that the integrand Xt in flH^ is given by 

xt = Eq [DtX I Tt] . (87) 

We now come back to the hedging equation flH^ . The fact that 6t = 6q 
is a solution to the homogeneous equation ( l85l) permits us to construct self- 
financed solutions of the in- homogeneous equation (l85l) . from solutions, which 
are not self-financed: 

Lemma 34 // 6 is an admissible portfolio (not necessarily self-financed) 
which satisfies [8^) . then there is a unique self-financing admissible portfolio 
9t such that the difference 9t — 9t is risk-free. It is given by: 

9t = atSo + 9t, (88) 

^t = -^AEQ[X\J't]-V,{9)]. (89) 

Proof. We here omit the argument u. Since the portfolio 9t — 9t is risk-free, 
it must have time to maturity 0, and the formula (188|) is true by definition. 
Substituting into equation ([BSD, and bearing in mind that (7^(0) = 0: 

{iPt(rt)'9ty =< 9t, ptai > = < at6o + 9t, Ptol > 
= at pt{0) (tI{0)+ < 9t, Ptai > 
= x\ Vz G I 

So 9t satisfies ( j85l) . It is then a hedging portfolio of X if Vt{9) = Eq[X \ JF^]. 
Substituting again (IHHl) and then ( 1891) . we get: 

Vt{9) = atVt{6o) + Vt{9) = atpt{0) + Vt{9) = Eq[X \ Tt\. 

If 9 is an admissible portfolio, then 9 is also admissible, since ||6'||p = ||^||p. 
■ 

By the lemma, the construction of a hedging portfolio for X is reduced 
to solve equation (!85l) in 9t{yj) for every (t, uS) , in such a way that 6^6?, i.e. 
9 is admissible. Any such solution 9 of this equation contains the risky part 
of the portfolio. 

To solve equation for given (t, cu), we have to know if XtioS) is in the 
range of the operator bti^u)'. The closure of the range of 6t(u;)' is equal to the 
orthogonal complement {)C{bt{uj)))^ of the kernel ICiptioj)) of bt{u}). 
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Consider the cases of I finite: The range TZ{{bt{w)y) is then closed, since 
it is finite dimensional. The kernel ]C{bt{uj)) is trivial iff the pt (uj) (t'I {uj) are 
linearly independent. So (64(0;))' is surjective and and there is a (non-unique) 
solution Otioj), for every Xt{ui), iff the pt iyj) o'l (^) linearly independent. 

Consider the cases of I infinite: The map {bt{uj))' from £'"'^([0, oo[ ) to 
is then never surjective. In fact, bt^uj) is a Hilbert-Schmidt operator, 
so it is compact. The adjoint is then also compact and since £^(I) is infinite 
dimensional, its range must be a proper subspace of This is the basic 

reason why there are always non-attainable contingent claims, when I is 
infinite. 

We have the following result (see Th.4.1 and Th.4.2 of [M] for the case 
I = N): 

Theorem 35 Let Dg = np>iLP(fi, P, J^). 

i) Ifl = N, then there exists X G Dg such that Vf{6) ^ X for all 6 G P^/. 

ii) Dq has a dense subspace of attainable contingent claims if and only if the 
operator cr^luj) has a trivial kernel a.e. {t,uj) G T x f2. 

Statement ii) says by definition that the bond market is approximately com- 
plete (notion introduced in [2] and [3]) if and only if (Xt^uj) has a trivial kernel 
a.e. 

In the sequel of this section, we are interested in the hedging problem 
for approximately complete markets, so we only consider the solution of 
the hedging equation ( l85l) in the case when cr^{uj) has a trivial kernel a.e. 
{t,Lu) G T X fi. 

Consider the case when I = N is an infinite and let i"^ = £^(I). To derive 
a condition under which flH31) has a solution and to derive a closed formula 
for one of the solutions, we rewrite the l.h.s. of ( l85l) using the notations 

It = CtPo, Bt{uj) = lt<Jt{uj) and r]t{uj) = S^'^ {pt{uj) / lt)9t{uj) . (90) 

Then 

{a^{u:))'pMOM = {a,{Lj)yk{pt{cu)/k)et{cu) = {ka,{u;)y{pt{Lu)/h)et{Lu) 
= {ka^{uj)rS-\pt{uj)/lt)et{uj) = iBtiuj)rvtiuj). 

The linear operator Bt{uj) is given, since po and a^{^JJ) are supposed given. 
Applying Theorem [2T] to the factor p//, it follows that equation ( l85ll is equiv- 
alent to find a progressive iJ^-valued process 77 satisfying the equation 

{Bt{uj)yrit{uo) = xt{uj), a.e. {t,uj) G T x fi. (91) 

We define the self-adjoint operator At{u!) in i"^ by 

At{uj) = {Bt{uo)yBt{uo). (92) 
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It is a fact of basic Hilbert space operator theory (cf. [T^) that the range 
n{{Btiuj))*) = n{{At{uj)y/'^). The solvabihty of each one of equations ([HSD 
and ( l9T|) is therefore equivalent to the existence of a progressive £^-valued 
process z satisfying 

{At{uj)y^^zt{uj) = xt{uj), a.e. {t,uj) eT xQ. (93) 

The kernel /C((A(o^))^/^) is trivial since /C((A(a^))i/2) = /C(A(cu)) = /C(5t(cu)) 
= {0}. Now, if xtico) E n{{Bt{uj)Y) then the unique solution of ^ is 
Zt{uj) = ({(At{uj)Y^'^)~^ Xt{uj) and a solution of (^1]) is given by 

r/iM = 5iH(AiH)"i/2xiH, (94) 

where St{uj), the closure of the operator Bt{uj){At{uj))~^^'^ , is isometric (cf. 
[IS]) from £2 to Let a be as in (fHSD and 

^ = a5o + ^ and = {lt/pt)Sr]t. (95) 

Then ^ is a hedging portfolio according to Lemma [HH 

In order to ensure that Xt{uj) of ( l85l) is in the range of (crjPt)(a;), we 
introduce spaces of vectors decreasing faster (for s > 0) than those of i"^. 
For s G M, let t''^ be the Hilbert space of real sequences endowed with the 
norm 

\\x\Usa = {J2{^ + ^yW\Y'- (96) 

Obviously f = f'^ and t''^ C £^'2, if s' > s. Although {At{iu))-^/^ is an 
unbounded operator in i"^ its restriction to ^'^'^ can be a bounded operator 
for some sufficient large s > 0, i.e. {At{uj))~^^'^i^''^ C This is the idea of 
our assumption, which will ensure hedgeability. However a precise formula- 
tion of this assumption must, as in the case of a finite of Bm., take care of 
integrability properties in (t,u). 

To consider also the case of a finite I, we define after obvious modifications 
the operator At^uj)) in £^(I) by formula (l92l) . In this case At{uj)) has obviously 
a bounded inverse. 

Condition 36 i) If Card{l) < oo, then there exists k G Dq, such that for 
allxE fil) : 

||x||^2 < k{uj)\\{At{uj)y/'^x\\i2 a.e. {t,uj) G T x fi. (97) 
a) IfI = N, then there exists s > and k G Dq, such that for all x G £^(I) : 
\\x\\i2 < k{uj)\\{At{uj)y/^x\Us,2 a.e. {t,u}) eT xQ. (98) 



39 



In the case of a finite number of Bm. Condition [32] i) leads to a complete 
market and one can choose a hedging portfolio such that it is continuous 
in the asset to hedge. To state the result let use introduce the notation 
Do(F) = r]p>iLP{n,P,J^,F), where F is a Banach space. Dq = Do(M). 

Theorem 37 (Finite number of random-sources, Card{T) < oo) 
If (0 (^f ConditionlSE is satisfied and if X e Dq, then the portfolio given by 
equation satisfies 6 G Pgf and Vf{6) = X. Moreover the linear mapping 
Dq3 X e P n Do(L2(T, E')), continuous. 

Proof. We only outline the proof of the theorem. Here i"^ = = is 
finite dimensional. 

Let X G Dq and let x be given by ( !83l) . First one proves (see Lemma 3.1 of 
[31]) that 

Do{F) = np>iLP{n,Q,J^,F). (99) 
Applying the BDG inequalities to equation flH^ it follows that 

xeD,iL\T,f)), (100) 

where x is progressively measurable. The definition of rj in and the 
condition fl^ give 

||?7t(t^)||£2 < kt{u;)\\xt{uj)\\e2. 

Inequality (ITOOl) then leads to r] e Dq{L^{T,E)). Using the definition ([95]) of 
6 we then obtain 

eeD,iL'iT,E')). (101) 

Since 6 satisfies equation fl85p by construction and since formulas fllOOp and 
fllOll) shows that 6 is admissible, the hypotheses of Lemma [3l] are satisfied, 
so 6* G Psf- This shows that is a hedging portfolio of X. 

All the linear maps X x ^ rj 6 are continuous in the above spaces, 
which also proves the claimed continuity of the map X 6. ■ 

The solution of the hedging problem, given by Theorem [371 is highly 
non-unique, since when Card{T) = m < oo then the kernel ]C{{af.pt){uj)) has 
infinite dimension. For instance there is a hedging portfolio consisting of 
m + 1 rollovers at any time. 

To state the result in the case of a infinite number of Bm., we first in- 
troduce spaces of contingent claims D^, smaller than Dq if s > and corre- 
sponding to that the integrand x in (l83l) takes values in More precisely, 
for s > let 

= {X G Dq I X G Do(L2(T,r'2)) where x is given by (^}. (102) 
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ConditionlHUlzz) leads to a D^-complete market, i.e. is a space of attainable 
contingent claims, is a dense subspace of Dq and is itself a complete 
topological vectorspace. This concept gives a natural frame-work to study 
existence and continuity of hedging portfolios. We have (see Theorem 4.3 of 

my- 

Theorem 38 (Infinite number of random-sources I = N) 

// (a) of ConditionlSE is satisfied and if X & D^, where s > is given by 
Condition then the portfolio given by equation l[9^) satisfies 9 G Pgf and 
Vf{e) = X. Moreover the linear map D^^X^OePn Dq{L'^{T, E')), is 
continuous. 

For the proof, which only uses elementary spectral properties of self-adjoint 
operators and compact operators, the reader is referred to 

A Malliavin-Clark-Ocone formalism was adapted recently in reference [6], 
for the construction of hedging portfolios in a Markovian context, with a 
Lipschitz continuous (in the bond price) volatility operator. This guaranties 
that the Malliavin derivative of the bond price is proportional to the volatility 
operator (formula (30) of [B]). Hedging is then achieved for a restricted class 
of claims, namely European claims being a Lipschitz continuous function in 
the price of the bond at maturity. 

References |H] and [27] studies the hedging problem in a weaker sense of 
approximate hedging, which in our context simply boils down to the well- 
known existence of the integrand x in the decomposition fl83p . 

6 Optimal portfolio management 

We now consider an investor, characterized by a von-Neumann-Morgenstern 
utility function [/, an initial wealth u, and a horizon T. The money is in- 
vested in a market portfolio, and the investor seeks to maximize the terminal 
(discounted) value Vf{9) of the portfolio. Transaction costs and taxes are 
neglected. The optimal portfolio problem is then to find an admissible self- 
financing portfolio 9 with Vo(^) = f , such that: 

r snvEp[U{Vf{e))] = Ep[U{Vf{e))] 

We will follow the now classical two-step approach (cf. [17], [28]) towards 
solving that problem. If the portfolio is self-financing and is worth v at time 
0, then, by the martingale property: 

Ep [^TVTm = V 
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where the random variable ^f, arising from Girsanov's theorem, was intro- 
duced earher in flTHj) . In general there can be several possible C,tj one for each 
7 satisfying the conditions of Definition [32l The first step (optimization) con- 
sists of finding for given 7, among jF^^-measurable random variables X such 
that Ep[^fX] = V, the one(s) that maximize expected utility Ep[U {X)]. 
This problem has in our setting a general solution X, given by Proposition 
Hni The second one (accessibility) consists in hedging one of the contin- 
gent claims X, obtained for the different 7, by a self-financing portfolio 6. 
This portfolio is then a solution of the optimal portfolio problem [Pq). By 
concavity, the final optimal wealth Vf{6) is unique. 



6.1 Optimization 

We consider, for a given 7 satisfying the conditions of Definition [221 the 
optimization problem: 

sup Ep [U (X)] 
Ep [^fX] = V 

X eL^ {n, Tf, p) 

We can rewrite it in a more geometric way, involving the scalar product in 
L2 (fi,^t,P): 

r sup U {X) dP 
J^^fXdP = i^f,X)L^=v 

( X e L^n,Tf,p) 

Problem (P) consists of maximizing a concave function on a closed linear 
subspace of L^. Assume there is a maximizer X. If the usual theory of 
Lagrange multipliers applies, there will be some A G M such that X actually 
optimizes the functional 



/ [U{X)-X^fX]dP 
Jn 



over all of L^. Maximizing pointwise under the integral, and bearing in mind 
that U is concave, we are led to the equation: 

U' (^X (uj)^ = X^f (uj) P-a.e., (103) 

which fully characterizes the solution X. Unfortunately this program cannot 
be carried through, for the function Ep [U {X)] has no point of continuity 
in unless U is bounded, so the constraint qualification conditions do not 
hold for problem (P), cf. [U]. We will therefore proceed by a roundabout 
way: use fll03p to define X, and then prove that X is optimal for a suitable 
choice of A. For this, we need some conditions on U. 
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Definition 39 The utility function U will be called admissible if it satisfies 
the following properties: 

1. f/ : M ^ {~oo} U M concave and upper semi- continuous 

2. there is some a G { — oo}U ] — oo,0], such that U (x) = —oo if x < a 
and U (x) > — oo if x > a 

3. U is twice differentiahle on the interval A = ]a, oo[; set B = U' [A) 
4- sup B = +oo; inf B = or inf B = — oo. 

5. U' : A ^ B is one-to-one, and there are some positive constants r, ci, C2 
and C3 such that its inverse I = [U']~'^ satisfies the estimate \I{y) \ + 
\yl'iy)\ < ci + C2 li/T + C3 \y\~'' for y e B. 

It follows from these assumptions that / is continuous and strictly decreasing, 
with: 

I (A) +00 when A inf S 
/ (A) a when A +00. 

We note that the estimate, in point 5) of Definition [39], is satisfied iff there 
exist C > such that 

\I{y)\ + \yI'{y)\<C {\yr+\yr), 
for all y E B. All usual utility functions are admissible: 
Example 40 

i) Quadratic utility; Set U (x) = jix — /i G M. Then a = —00, and 
U' [x) = yU — X, so that B = M. and I {y) = fi — y. The estimate is satisfied 
with r = 1. 

a) Exponential utility; Set f/ (x) = 1 — ^ exp {—fix) , /i > 0. Then a = —00, 
and U' (x) = exp (— /ix) , so that B =]0, oo[ and I (y) = —j^\n{y). The 
estimate is satisfied for any r > 0. 

Hi) Power utility; Set U (x) = ^x^ for some fi < 1 and (note that fi 
may be negative). Then a = 0, and U' (x) = x^~^ , so that B =]0, 00 [ and 
I (y) = y^/^t^~^\ The estimate is satisfied with r = 

iv) Logarithmic utility; Set U (x) = Inx. Then a = and U' (x) = ^, so that 
B =]0, 00 [ and I (y) = -. The estimate is satisfied with r = 1. 
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Take some \ E B and a 7 satisfying the conditions of Definition [321 and 
define a random variable Xx by: 

Xx is jFj^-measurable. In addition, we have: 
Lemma 41 G (^2, JF^^, P) for every p > 1. 

Proof. Since U is admissible, we know from condition 4 that, for some r > 
we have: 

\nx^f)r<{ci+c2\x^T\''+c3\x^frY 
< A;i + A;2iAriefr + ^3iAr^^ leti"''' 

and the right-hand side is integrable, for we know that G {Q,J^f,P) 
for every s G M. ■ 

Lemma 42 Let v E A. There is a unique \ E B such that Ep [X-^^f \ = v 

Proof. Consider the map (p : B ^ M. defined by Lp (A) = Ep [Xx^f] = 
Ep [I (X^f) ^f]. Since > P-a.e., and / is strictly decreasing, ip is strictly 
decreasing. Using the Lebesgue dominated convergence theorem, we find 
that it is continuous. Using Fatou's lemma, we find that: 

• (fi (A) — ^ +00 when A — > inf P 

• limsup (p (A) < a when A +00 

Since v E A, it follows that there is a unique A such that (p (jCj = v. m 

Denote X-^ by X. We now conclude: 

Proposition 43 X is the unique solution of problem (P). 

Proof. Let us show that X is indeed a solution of problem (P). Uniqueness 
follows from the strict concavity of U. 

We have shown that X is in L^, and Ep[X^f] = v, so X satisfies the 
constraints. Take another X G such that Ep[X^f] = v. Since U is 
concave, we have: 

U (X {00)) <U(^X (cu)) + (X {00) - X {uj))U' (x (cu)) P-a.e. 



44 



By definition, U' yX {uj)j = XC^f{uj). Substituting into the inequality and 
integrating, we get: 

[ U{X)dP< [ u(x)dP + X I {X- X)ifdP 
Jn Jn ^ Jn 

and the last term vanishes because it is just X{v — v). So X is indeed an 
optimizer, and the result follows. ■ 



6.2 Hedging 

Once the solution X of the optimization problem (P) is found, for a given 
7, the question is whether it can be hedged by a self-financing portfolio 6, 
so that Vf{9) = X. We note that, if there exists such 9 G P^/, then it is a 
solution of (Po)- In fact, let 6 G P^/ and Vo{6) = v and set X = Vf{6). It 
follows from (P) that 



Ep [U{Vf{9))] = Ep [U{X)] < Ep U{X) 



Er 



UiVfiO)) 



so 6^ is a solution of (Po). 



6.2.1 Deterministic case 

In this paragraph, we shall use the general hedging results of ^to solve this 
problem, in the case when the m and a, are deterministic (i.e. they do not 
depend on u). 

Under these conditions, there can be several 7 that satisfy the conditions 
of Definition [HI] and some 7 can even be non-deterministic. However, as we 
have supposed that the market is strongly arbitrage free, so equation (17^ 
has a solution, we can choose 7 to be the unique solution with the property of 
being orthogonal in i"^ to the kernel of the volatility operator. More precisely, 
we choose the unique 7 such that 

(7^,x)£2 = 0, Wxef{I) s.t. atx = 0. (104) 

The 7 defined by this condition is deterministic. In the sequel of this para- 
graph 7 is given by (11041) . In that case, it follows from formula ( 1761) that 
is Malliavin differentiable. It follows from formula (l86ll that the partial 
derivative with respect to is given by: 
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and X = I {\^f^ is Malliavin differentiable as well, with: 
The Clarke-Ocone formula now reads: 



X = I J-q] + V / Eq [D^^tX I dWl (105) 

^ - ^ E [^T/'lACf ) I -^t] rfW^/ (106) 



We then write the equation ( l85|) for the hedging portfolio ^, and we substitute 



the Clark-Ocone formula for x\ {uj): 

h{u)'et{uj) = -xeq [ifi'ix^t) I J^t] If (107) 

This equation has a solution iff 7^ is in the range of ht{uj)' . Since a is determin- 
istic, this condition simplifies. In fact, let It and Bt be given by ( l9Qll . which 
here both are deterministic, and let q{t,uj) = pt{uj)/lt. Then the expression 
(IH2D of bt{iu)' give: 

(6iH'^i(a;))^ =< ^i(a;) , >=< 9t{u)q{t,u) , kal >= {B'Jt{u))\ 

where ftil^) G E~'' is given by ftif^) = q{t,uj)6t{uj). So, equation f ll07p has a 
solution iff 7^ is in the range of B[. This is always true when I is finite, since 
then the range of B'^ is equal to the orthogonal complement of the kernel of 
(jf (we remember that pti'^^x) > for x > 0). When I = N, then the range 
is only a strictly smaller dense subset. 
We are lead to following condition 

Definition 44 We shall say that the market satisfies condition (C) if there 
exists a deterministic portfolio 6^ which is admissible and satisfies B[6^ = 74, 
i.e. 

<e^,{CtPo)ai> =Yt, (108) 

for each i El and t. 

Condition C is then equivalent to 7^ G 1Z{B[), the range of B[. In the case 
when I is finite, there is never uniqueness in the choice of 9^. 

In the case when I is finite, we know that condition (C) is satisfied and 
it can easily be verified, with n elements say, by picking n maturities < 
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< ... < 5^ and by seeking 9^ as a linear combination of rollovers: 9^ = 

xl6si- Condition fllUSI) then reduces to a system of n linear equations with 
n unknowns which determines the x]. 

In the case when I = N, condition (C) may not be satisfied. We will be 
content with reminding that the left-hand side of equation f llOSp is meaning- 
ful, since (CtPo) al belongs to the space E'^ . 

If condition (C) is satisfied, equation (I107p becomes: 

< 9,, p,ai > = -XEQi^fl'Cx^f) I J-,] < 91 > 

and an obvious solution 9t = 9t (the risky part of the optimal portfolio) is 
given by: 

9, = -XEgl^fl'Cx^f) I J^t] ^ 9^ 

Applying Lemma 1341 with x defined by ( l89l) . we obtain a hedging portfolio 
9 = x6q + ^ of X, where 9 is as above, and: 



Xt 



Pt (0) 



(^EQ[l(x^f^ \J^t]-<9t,pt>). 



To sum up, in the case when the and the al,i G I, are deterministic, 
with al (0) = 0, with condition (C) and equation fl74p satisfied, an optimal 
admissible and self-financing portfolio is given by 

9t = XtSo + 9t, where 9^ = yt 9^ (109) 

Pt 

and where the coefficients Xt and yt are real-valued progressively measurable 
processes given by 

yt = -EQ[X^fI'CX^f)\J't] (110) 

= ipmr' (EQinkf) i -^d -yt< o'^ , >) . (m) 

This leads immediately to a mutual fund theorem: whatever the utility 
function and the initial wealth, the optimal portfolio at time t is a linear 
combination of the current account So and the portfolio / t— 9^, -^^/ >, 

i.e. the portfolio ^^9^. This portfolio is in general not self-financed, so it 
can not be given the status of a market portfolio. However we can easily 
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reformulate the result with a self-financed portfolio. In fact, chose an ad- 
missible utility function, with a = 0, according to Definition [23 For this 
utility function, let 9 be the optimal portfolio given by (11091) . with unit ini- 
tial wealth. Obviously is a linear combination of 5q and 6^. This gives 
us: 

Theorem 45 (Mutual fund theorem) The optimal portfolio G has the 
following properties: 

i) Q is an admissible self-financing portfolio, with unit initial value, i.e. 
< ©0 ; >= 1; (in-d the value at each time t E T is strictly positive, i.e. 
<Qt,Pt » 0. 

a) For each utility function U, admissible according to DefinitionlEB o-nd each 
initial wealth v G ]a, oo[, there exist two real valued processes c and d such 
that if 9t = CtSo + dtQf, then 9 is an optimal self financing portfolio for U, 
i.e. a solution of problem (Pq). 

6.2.2 Stochastic m and a 

We shall here concentrate on the case of an approximately complete market, 
which is equivalent to that the volatility operator is non-degenerated. In 
fact, according to Hi) of Theorem [SHI the market is approximately complete 
if and only if (J^{uj) has a trivial kernel a.e. {t,uj) G T x fi. We remind that 
the market of price process 7 is unique in this case. 

In the case of a finite number of Bm. we obtain easily from Lemma |4T] 
and Theorem [S7I the following result (see Theorem 3.6 of [TU]): 

Theorem 46 Let I be a finite set, let U be admissible in the sens of Defini- 
tion{3M and let i) of Condition{3E be satisfied. The problem (Pq) then has a 
solution 9. One solution 9 = ado + 9 E Pgf is given by ( [P3]) . 

In the case of an infinite number of Bm. we shall impose Malliavin differen- 
tiability properties on the market price of risk 7. To this end we introduce 
the space Dj, for s > by 

Dl = {X eDqIDX e Do(L2(T,r'2))}. (112) 

We can now state a result in the case of an infinite number of Bm., quite 
analog to the case of a finite number of Bm. (see Theorem 4.5 of [M]): 

Theorem 47 Let I = N, let U be admissible in the sens of Definition let 

ii) of Condition [3^^ be satisfied and let ln(^2^) G D^, where s > Q is given by 
ii) of Condition\3^ The problem (Pq) then has a solution 9. One solution 
9 = aSo + 9 E Psf is given by f95\) . 
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Proof. We only consider the case of U' > 0, since the case of U'{x) = for 
some X is so similar. Let the hypotheses of the theorem be satisfied. The 
portfolio ^ is a solution of equation (Pq), if 9 E Psf and if it hedges X given 
by Proposition HHl (See Corollary 3.4 of [in])- It is enough to verify that 
Theorem [381 applies to X = I{XC,f) for a certain given A > 0. 

/ is C\ so DtX = X^fip'{X^f)Dtln{^f). Since In(^f) e D], this gives 
||M||i2(x,£s,2) = \X^f^'iX^f)\ ||^ln(^f)llL2(T/=.2)- The inequality in 5) of Def- 
inition ESI gives ||M|U2(T,,.,2) < CiiX^fY + iX^f)-P)\\D\n{^f)\\L2^j^,s,2), for 
some p>l. Condition ^ of Definition [32] shows that (A^^)^ + (^-^f)"^ ^ 
L^(fi,P), for all g > 1. By hypothesis ||D ln(^'f)||£,2(T,£s,2) G Dq, so Holder's 
inequality now gives that ||L'X||i2('ir,^s,2) e Dg, i.e. DX G Dg(L^(T, By 
Lemma im X G Dg. It follows that X G D]. We can now apply Theorem [38l 
which proves the existence of 9. ■ 



6.2.3 Examples. 

We now give some examples of optimal bond portfolios for logarithmic and 
quadratic utility functions U. Other examples can be found in [TU] . 

First we assume the drift function rrit and the volatility operator at to be 
deterministic. We shall therefore suppose that the market satisfy condition 
(C), of Definition HU so the market prices of risk 7 is deterministic and 
satisfy condition fllUSp . We shall derive the optimal portfolio directly, going 
through the steps leading to the general solution fll09p . 

Secondly we study the general case of stochastic drift function rut and 
volatility operator at for the logarithmic utility function. 

The final optimal discounted wealth is X = I{X^f). The corresponding 
optimal discounted wealth process Fis given by Yt = Eqll^X^f) | JF^]. The 
initial wealth Yq = v determines A by the equation 

v = Yo = EQ[lCX^f)]. (113) 

We recall that {pt)~^CtPo G E'^ a.s and that pt{0) > a.s. 

Logarithmic utility (deterministic m and a) Let 

U{x) = ln(x). (114) 
We have I{x) = 1/x, and X = (A.^f)~\ so that equation (11131) gives: 
V = EQ[1/Cx^t)] = Ep[^f/Cx^f)] = 1/A. 
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Then using the expression flTUj) for C,t and we have: 

} = exp(-i fY.Hfds+ tj2^ldW:]. (115) 

The right-hand side is a Q-martingale, then so is 1/C,t- It follows that the 
optimal discounted wealth at t is 

Y, = EQ[lCX^f)\:Ft] = ^ = y. 
Since rf(l/6) = T.ieii^t/^t)dWi and X = Yf, it then follows that: 

The hedging equation ( 185|) and the above formula give: 

Vz G I, < (u) , Pi (.7) (cu) > = -^7l (117) 

By condition (C) we find a portfolio 6''^ satisfying 7^ = < 6'° , (CtPo) al >, so 

Y, = <{C,po)elal>. (118) 
Substituting this expression of 7 into (11171) we obtain: 

Vz G I, < Pi (^) (uj) - -f-(Apo) {00) > = 0. (119) 

6(^) 

One solution of this equation is obviously given by ^ = ^, where 

Otiuj) = miuj) e^, ytiuj) = -f-. (120) 
Pt{^) 6(^) 

The discounted value of 9 at time t in state u is then 

mm^) =< e),,p, >= -f- < elcpo > . (121) 

The optimal portfolio 9 is now obtained by using Lemma 1341 9t = XtSo + 9t, 
where 

Pi(0) 6 

As it should, the discounted value of 9 is then Vt{9) = Yt = v/^f 

We note the following useful property: the ratio of the investment in 
bonds with time to maturity 5 > to the total investment is deterministic. 
In fact this ratio is simply price at t = 0, of a zero-coupon bond with time 
to maturity S + t : 

9t{S,uj) pt{S,uj) 



Po{S + t). (123) 
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Quadratic utility (deterministic m and a) Let the utility function be: 

U {x) = /ix — -x^ 
As in i) of Example HQl we find that 

Hy) = f^-y- 

The final discounted optimal wealth is X = I{X^f), so 

X = /i — XC,f. 
We determine A by the condition: 



v = E, 



Q 



X 



Er 



/i - XEq [^f] 



Set 



\ ^Jo Jo 
Then Z is a martingale with respect to Q and formula ( 1771) gives 

We have, by substitution into (11241) : 

V = /i - XEq [^f] = /i - A exp / y^illfd. 

This gives 

X = {^i-v)exp(^- £^^{jlfds^. 

It now follows from (I125p that 

X = fi — X^f = fi + [v — fi) Zf 
and the optimal discounted wealth at t is 

Yt = EQ[lCX^f)\J^t] = f^+{v-fi) Z,. 



(124) 



(125) 



(126) 



(127) 
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Since dZt = —ZtYlinzilldWl, we have that 

so the hedging equation reads (see (l85ll ): 

Vz G I, < Ot {uo) , pt (^) {uo)>= - (/i - Yt{uo)) 7!- (128) 

As usually, condition (C) gives a portfolio 6^ satisfying 7^* = < 6^° , (CtPo) >, 
which together with (11281) gives: 

Vz G I, < Pi (u) Ot [u] + (Fi(c^) - /i) (Apo) ^1 (^) > = 0. 

One solution of this equation is 6' = ^, where 

U^) = y,{co) 9l y,{u) = /i - Y,{uj). 

9 gives the risky part of the optimal portfolio. 

Applying Lemma IMl we obtain the optimal portfolio 9t = XtSo + 9t, where 

= iptm~\yit) - (/i - Y{t)) < ^° , CtPo >). (129) 



Logarithmic utility (stochastic m and a) We assume that the condi- 
tions of Definition [32] are satisfied. We chose 7^(0;) to be orthogonal to the 
kernel of (t^{uj), a.e. {t,uj). This 7 satisfies the conditions of Definition [321 
Formulas (I114l) - (lll7p then still hold true. As in the discussion preceding the 
condition (C), of Definition [H] it follows that 7t(cu) is a.s. in the closure of 
the range of Bl{uj). Therefore, in this example, the natural generalization of 
the condition (C) to the stochastic case is simply to impose the same condi- 
tion (I108p of Definition [m to be satisfied with a stochastic portfolio 9^ G P. 
Formulas (I118p - (I123I) are then also true statements and it follows using The- 
orem [21] that 9 G Pgf- In particular the ratio of the investment in bonds with 
time to maturity S* > to the total investment is deterministic. 



6.3 The H-J-B approach 

When rrit and al are given functions mt{pt) and (Jl{pt) of the price pt-, for 
every t, then the optimal portfolio problem (Pq) can be considered within 
a Hamilton- Jacob i-Bellman approach. In this subsection we illustrate this 
approach, without being rigorous and we suppose that the utility function U 
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satisfies the conditions of Definition [311 For notational simplicity we exclude 
the price argument in rrit and a^. 

The optimal value function, here denoted by F, then only depends of time 
t, of the value of the discounted wealth w and the discounted price function 
f & oi Zero-Coupons at time t : 

F{t,wj) = sn^{E[u{yf{e)) I vm = w,vt = f]\ee p,/}. 

The derivative DG{f] g) of a function 3 f G{f) in the direction g E E"^, 
is as usually defined by 

Suppose that G is C\ Writing DG{f) for the map g ^ DG{f; g) and D'^GW) 
for the map gi x g2 ^ DG{f\ gi, g2), we have that DG{f) is a linear contin- 
uous form on E^ and D'^G{f) is a bi-linear continuous form. 

Let us first consider the case of a volatility operator a with trivial kernel, 
i.e. for every strictly positive price (function) f E E^, the kernel of the linear 
map (Tj : — > E'^ is trivial a.s. According to Definition [32] there is then a 
unique market of price process 7. Define the Hamiltonian H{t,w, f,x) by: 

H(t, w, f,x) = J2 ^'{t^ fHg^ (t^ f) + DF{t, w, f;df + Yl <^lf) 



(130) 



In that formula, x = (a;*)^gj G £^ is the control, which is related to the optimal 
terminal wealth by formula (l83i) . A control x is called admissible if 

x\t,Vtie),pt)=<et,ptcTi> (131) 

for all 6 G Ps/. In other words, x* can be interpreted as the value invested 
in the z-th source of noise. Using the Ito formula, one derives the (formal) 
HJB equation: 

dF 

— {t, w, f) + sup H{t, w, /, x) = 0, (132) 
with the boundary condition 

F{f,w) = U{w). (133) 
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The optimal control £, solution of the optimization problem 

sup H{t,w J, x), 

X 

is given by 

Now, substitution of H(t,w, f,x(t,w, f)) into equation f ll32p gives: 



(t, w, f) f — (t, w, f) + DF{t, w, f; df + mj) 



dw'^ ' ' \ dt 



dF dF 



(135) 



Once the solution F of (11351) . with boundary condition (I133p . is found, the 
optimal control x is given by (11341) . Any optimal portfolio 6 is then a solution 
of the equation: 

x\t, Vti9),pt) =< Ot , ptai >, V z G I, t G T. 

Next we consider the case of a volatility operator, which does not necessar- 
ily have a trivial kernel. Once more we define the Hamiltonian i?(t, w, /, x, 7) 
by formula (I130p . which now also depends on the control 7, a £^(I) valued 
function of {t,w,f). A control (x,7) is admissible if condition (I13ip is sat- 
isfied and if the conditions of Definition [32] are satisfied, so writing out the 
price argument / G -E* in rrit and al : 

mtif) = atifMwJ). (136) 

The optimal control 7 is determined by conditions (I13ip and (I136p . This 
can be seen as follows. Let 7"'"(/) be the unique solution of (11360 such that 
j^if) is in the orthogonal complement {)C{at{f)))^ of the kernel /C((Tt(/)), 
let d = 7 — 7-*- and let Pt{f) be the orthogonal projection on IC{at{f)). 
Condition (113 ip implies that x G {IC{at{f)))^ ■ According to (I134p . this can 
only be satisfied if 

dF 

7 = 7-L + a and at{wJ)— = Pt{f)vt{wJ), (137) 

where i'l{w, /) = {D^){t, w, f; crlf). So in the general the case of a volatility 
operator, which does not necessarily have a trivial kernel, the H- J-B approach 
leads to the equation (I135p . with 7 replaced by 7 defined by formula (I137p . 
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In the case when rrit and al are independent of pt, then the f * are inde- 
pendent of /, 7 = 7-*- and the above equations simphfy: 



dF 92 F 



dt dw^ 

with the boundary condition 

F{T,w) = U{w), w eR. 
Each self financing portfolio 6 G Ps/, such that 



\ dw J \ dw'^ 

where w = Vt{6), is then a solution of problem (Pq)- The solutions in the 
examples in §6.2.31 as well as the general solution fllOOp for deterministic m 
and a, are easily obtained by solving these equations. 



A Appendix 

In this appendix, we reproduce results (proved in the appendix of [ID]), used 
in this article, concerning existence of solutions of some SDE's and esti- 
mates of these solutions. The notations T = [0, T], W^, I and (f2, P, JF, A) are 
defined in §3.11 Through the appendix m and cr*, i G I, are ^-progressively 
measurable £''*-valued processes satisfying 



/ (||mt||£;s ||crt|||s)(it < oo,a.s. (138) 
Jo 



'0 iei 
The i^'^-valued semi-martingale L is given by 



L{t)= I [msds + ^aidWi), if < t < T (139) 

and by L[t) = L{T), if t > T. We introduce, for t > 0, the random variable 
Kt)=t+ [ {\\m,\\Es + J2\\<\\l^)ds, ifO<t<f (140) 



55 



and fi(t) = t — T + /i(T) if t > T. /i is a.s. strictly increasing, absolutely 
continuous and on-to [0, oo[. The inverse r of /i also have these properties 
and r(t) < t. For a continuous iJ'^-valued processes Y on [0, T] we introduce 

p,{Y) = {E[snp \\Y{r{smUy/\ (141) 
se[o,t] 



for t G [0, oo[, where we have defined Y{t) for t > T by Y{t) = Y(T). We 

'^se[o,t] 



note that pt{Y) < {E[sup,^,Q^t] since r(t) < t. 



Lemma 48 // condition I^138\) is satisfied and if Y is an A-progressively 
measurable E'^ -valued continuous process on [0,T], satisfying pt{Y) < oo, for 
all t > 0, then the equation 

X{t) = Y{t) + [ Ct-sX{s){msds + J2 '^id^D^ (142) 

t G [0,T], has a unique solution X, in the set of A-progressively measurable 
E"^ -valued continuous process satisfying: 

[ i\\Xit)\\Es + \\Xit)mt\\Es + V \\Xit)al\\%)dt < oo a.s. (143) 
Moreover this solution satisfies: 

i) If /o^(||?^f lls^+i + X]jei 11*^1 II B'i+i)'^^ < andY is a continuous E^^^-valued 
process with pt{dY) < oo, for all t > 0, then X is a continuous E^^^ -valued 
process. 

ii) If (z) is satisfied and if Y is a semi-martingale, then X is a semi- 
martingale. 

Hi) IfY is -valued, then X is H'^ -valued. 



The next lemma establish conditions under which the solution of equation 
(11421) is in L^, p G [0, oo[ . The notation £ was introduced in 

Lemma 49 Let condition U38\) be satisfied and let (i) 

E[exp{p / {\\mt\\Es + \Wl\\ls)dt)] < oo, 
"^0 iei 

for each p E [1, oo[. Suppose that Y in Lemma\48\ satisfies (ii) 

E[sup\\Yit)\\U<oo, 
teT 
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for each p E [1, oo[. Then the unique solution X of equation (TJ^ in Lemma 
^ satisfies 

E[sup\\X{t)\\%] < oo, Vpe [l,oo[. (144) 

Moreover if (Hi) 

E[{ [ {\\mt\\Es+. + J2 \Wl\\%+^)dm < oo 
"^0 m 



and (iv) 



E[snp\\Y{t)\\%+,] < oo, 



for each p G [1, C)o[, then also 

E[sup||X(t)||^.+,] < oo, Vpe [l,oo[. (145) 

In particular, estimates il44\) '^^^ ( [i^5| ) applies to X = S{L). 
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Note added in the proofs: Since tlie preparation of this paper, the optimal 
bond portfolio problem has further been studied in various directions: 

1 The reference De Donno, M. and Pratelli, M.: A theory of stochastic 
integration for bond markets, Ann. Appl. Probab. 15, 2773-2791 
(2005) considers the optimal bond portfolio problem in a more general 
semi martingale bond market. Existence of optimal wealth strategies 
is established and existence of optimal portfolios is studied. 

2 The reference Ringer, N. and Tehranchi, M.: Optimal portfolio choice 
in the bond market, Finance Stoch. 10, 553-573 (2006) considers 
the optimal bond portfolio problem in a Markovien setting of local 
volatility operators with full range and which are globally Lipschitzien. 
More precisely it is assumed, with our notations and limiting us to the 
time homogeneous case, that the function C : E ^ T-lS{i'^,E), where 
C(/) = f<^{f)i is globally Lipschitzien and that for all strictly posi- 
tive f E E the closure of the range lZ{C{f)) is the subset of elements 
g E E such that 5^(0) = 0. If moreover (the unique) market price of 
risk is globally Lipschitzien then they establish the existence of a solu- 
tion to the optimal portfolio problem. We note that the proof of this 
boils down to the verification of properties of the Malliavin derivative 
of ln(^f ) as was already the case in Theorem 4.5 of [M] (see Theorem 
|47|) . We also note that their Gaussian example, of course satisfies our 
condition (C) of Definition HU so it is covered by our treatment. 
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